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We  propose  a  method  that  enables  strong,  coherent  coupling  between  individual  optical  emitters  and  elec¬ 
tromagnetic  excitations  in  conducting  nanostructures.  The  excitations  are  optical  plasmons  that  can  be  local¬ 
ized  to  subwavelength  dimensions.  Under  realistic  conditions,  the  tight  confinement  causes  optical  emission  to 
be  almost  entirely  directed  into  the  propagating  plasmon  modes  via  a  mechanism  analogous  to  cavity  quantum 
electrodynamics.  We  first  illustrate  this  result  for  the  case  of  a  nano  wire,  before  considering  the  optimized 
geometry  of  a  nanotip.  We  describe  an  application  of  this  technique  involving  efficient  single-photon  genera¬ 
tion  on  demand,  in  which  the  plasmons  are  efficiently  outcoupled  to  a  dielectric  waveguide.  Finally,  we  analyze 
the  effects  of  increased  scattering  due  to  surface  roughness  on  these  nanostructures. 
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I.  INTRODUCTION 

In  recent  years,  there  has  been  substantial  interest  in 
nanoscale  optical  devices  based  on  the  electromagnetic  sur¬ 
face  modes  (surface  plasmons)  associated  with  subwave¬ 
length  metallic  systems.  Surface  plasmons1  are  electromag¬ 
netic  excitations  associated  with  charge  density  waves  on  the 
surface  of  a  conducting  object.  The  unique  properties  of  plas¬ 
mons  on  nanoscale  metallic  systems  have  produced  a  num¬ 
ber  of  dramatic  observed  effects,  such  as  single  molecule 
detection  with  surface-enhanced  Raman  scattering,2,3  en¬ 
hanced  transmission  through  subwavelength  apertures,4,5  and 
enhanced  photoluminescence  from  quantum  wells.6  There  is 
also  considerable  interest  in  these  systems  in  applications 
such  as  biosensing,7  subwavelength  imaging,8,9  and 
waveguiding  and  switching  devices  below  the  diffraction 
limit.10-13  Such  subwavelength  waveguiding  of  plasmons  in 
metallic  nanowires  has  been  observed  in  a  number  of  recent 
experiments.14-16 

At  the  same  time,  spurred  in  part  by  rapid  developments 
in  the  fields  of  quantum  computation  and  quantum  informa¬ 
tion  science,  there  has  been  strong  interest  in  exploring  new 
physical  mechanisms  that  enable  coherent  coupling  between 
individual  quantum  systems  and  photon  fields.  Such  a 
mechanism  would  enable  quantum  information  to  be  passed 
over  long  distances  and  long-range  interactions  between  sys¬ 
tems.  These  features  are  not  only  essential  for  quantum 
communication17,18  but  would  also  facilitate  the  scalability  of 
quantum  computers.19  The  required  coupling  between  emit¬ 
ters  and  photons  is  difficult  but  has  been  achieved  in  a  num¬ 
ber  of  systems  that  reach  the  so-called  strong-coupling 
regime  of  cavity  quantum  electrodynamics  (QED)  20-22  Re¬ 
cently,  several  approaches  to  reach  this  regime  on  a  chip  at 
microwave  frequencies  have  been  suggested23-25  and  experi¬ 
mentally  observed,22  which  utilize  coupling  between  emitters 
and  modes  of  superconducting  transmission  lines.  A  key  fea¬ 
ture  of  these  transmission  lines  is  the  reduction  of  the  effec¬ 
tive  mode  volume  Veff  for  the  photons,  which  in  turn  results 
in  a  substantial  increase  of  the  emitter-field  coupling  constant 


grc  1/Vveff-  Realization  of  analogous  techniques  with  optical 
photons  would  open  the  door  to  many  potential  applications 
in  quantum  information  science,  and  in  addition,  lead  to 
smaller  mode  volumes  and  faster  interaction  times. 

In  this  paper,  we  describe  a  method  that  enables  strong, 
coherent  coupling  between  individual  emitters  and  electro¬ 
magnetic  excitations  in  conducting  nanostructures  on  a  chip 
at  optical  frequencies,  via  excitation  of  guided  optical  plas¬ 
mons  localized  to  nanoscale  dimensions.  The  strong  coupling 
occurs  due  to  the  subwavelength  confinement  and  small 
mode  volumes  associated  with  the  surface  plasmon  modes. 
We  show  that  under  realistic  conditions  optical  emission  can 
be  almost  entirely  directed  into  these  modes,  in  a  manner 
similar  to  the  Purcell  effect  in  cavity  QED.  We  first  examine 
the  case  of  a  cylindrical  nanowire,  a  simple  geometry  where 
the  relevant  physics  can  be  understood  analytically,  before 
considering  the  more  optimal  geometry  of  a  conducting  nan¬ 
otip.  Defining  a  Purcell  factor  P  =  rpl/rother  characterizing 
the  ratio  of  the  spontaneous  emission  rate  into  the  plasmon 
modes  over  emission  into  other  channels,  we  show  that  val¬ 
ues  exceeding  P~  103  are  possible  in  these  systems,  limited 
only  by  metal  losses  at  optical  frequencies.  Because  of  these 
losses,  the  plasmon  modes  themselves  are  not  suitable  as 
carriers  of  information  over  long  distances.  However,  we 
show  that  the  plasmon  excitation  can  be  efficiently  converted 
into  a  photon  in  a  nearby,  evanescently  coupled  dielectric 
waveguide,  as  illustrated  schematically  in  Fig.  1.  This  can  be 
used,  e.g,  to  create  an  efficient  single-photon  source  or  as 
part  of  an  architecture  to  perform  controlled  interactions  be¬ 
tween  distant  qubits.  The  achievable  coupling  between  the 
plasmon  and  waveguide  systems  can  be  much  stronger  than 
the  plasmon  dissipation  rates,  and  we  find  that  single-photon 
generation  efficiencies  exceeding  95%  are  possible  for  the 
simple  geometries  considered  here. 

This  paper  is  organized  as  follows.  In  Sec.  II,  we  calculate 
the  mode  structure  of  a  conducting  nanowire  surrounded  by 
some  positive  dielectric  medium.  We  show  that  the  nanowire 
supports  one  fundamental  plasmon  mode  with  significantly 
reduced  phase  velocity,  which  is  tightly  localized  on  a  scale 
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FIG.  1.  (Color  online)  (a)  An  emitter  coupled  to  a  nanowire  is 
optically  excited  and  decays  with  high  probability  into  the  plasmon 
modes  of  the  nanowire.  A  single-photon  source  is  created  by  eva- 
nescently  coupling  the  nano  wire  to  a  nearby  dielectric  waveguide 
over  a  length  Lex.  The  single-photon  source  can  potentially  be  uni¬ 
directional,  e.g.,  by  capping  one  end  of  the  waveguide  with  a  re¬ 
flective  surface,  (b)  An  internal-level  scheme  that  allows  for  shap¬ 
ing  of  the  outgoing  single-photon  pulses.  An  emitter  that  starts  in 
state  | s)  is  coupled  to  excited  state  \e)  via  a  time-dependent  external 
control  field  £l(t).  We  assume  that  the  excited  state  \e)  is  coupled  to 
state  | g)  via  the  plasmon  modes,  causing  \e)  to  decay  into  | g)  with 
high  probability,  while  simultaneously  generating  a  single  photon  in 
the  plasmon  modes.  The  shape  of  the  photon  wave  packet  is  deter¬ 
mined  by  £l(t).  (c)  A  similar  scheme  for  single-photon  generation 
using  an  emitter  coupled  to  a  nanotip  instead  of  a  nano  wire.  Note 
that  this  scheme  is  naturally  unidirectional,  as  the  generated  plas- 
mons  propagate  in  a  single  direction. 

~R  around  the  wire  surface.  We  also  calculate  the  dissipa¬ 
tion  rate  of  the  fundamental  mode  as  it  propagates  along  the 
nano  wire,  due  to  metallic  losses.  In  Sec.  Ill,  we  calculate  the 
emission  properties  of  a  dipole  emitter  near  the  nanowire  as 
a  function  of  emitter  position  and  wire  radius.  We  show  that 
under  certain  circumstances,  emission  into  the  guided  plas¬ 
mon  modes  is  greatly  enhanced  over  decay  into  radiative  and 
nonradiative  channels.  In  fact,  when  optimized,  the  probabil¬ 
ity  of  emission  into  the  plasmon  mode  approaches  almost 
unity  for  small  R  and  is  limited  only  by  inherent  losses  in  the 
metal.  While  the  nanowire  is  a  simple  system  to  understand, 
we  note  that  there  is  an  inherent  trade-off  present  in  this 
geometry,  as  the  small  wire  sizes  that  yield  strong  coupling 
also  result  in  more  rapid  dissipation  of  the  plasmons  as  they 
propagate.  Therefore,  in  Sec.  IV,  we  consider  a  better- 
optimized  system,  a  conducting  nanotip.  We  show  that  the 
nanotip  can  significantly  reduce  the  effects  of  propagative 
losses  even  while  preserving  large  coupling  strengths.  In  Sec. 
V,  we  consider  the  problem  of  outcoupling  the  plasmon 
modes,  and  study  in  detail  the  interaction  between  the  plas¬ 
mon  modes  of  our  nanostructures  and  the  guided  modes  of  a 
nearby  dielectric  waveguide.  We  show  that  the  plasmon 
modes  can  be  efficiently  outcoupled  to  the  waveguide,  and 
we  propose  an  architecture  for  efficient  single-photon  gen¬ 
eration  on  demand  based  on  a  tiered  emitter/nanostructure/ 
waveguide  system.  We  calculate  the  expected  efficiencies  for 
single-photon  generation,  taking  fully  into  account  the 
propagative  losses  of  the  plasmons,  the  finite  Purcell  factors 
governing  the  interactions  with  the  dipole  emitter,  and  the 
nonunity  coupling  efficiency  between  the  plasmon  and  wave¬ 
guide  modes.  In  Sec.  VI,  we  consider  the  effects  of  possible 


imperfections  to  the  system,  in  particular,  the  adverse  effect 
of  surface  roughness  on  the  conducting  nanostructures.  In 
general,  surface  roughness  can  lead  to  radiative  scattering  of 
plasmons  as  well  as  increased  nonradiative  dissipation, 
which  results  in  larger  losses  as  the  plasmons  propagate 
along  the  structure.  We  calculate  the  effects  of  these  two 
processes  and  find  only  moderate  increases  in  the  total  loss 
under  reasonable  parameters.  Finally,  in  Sec.  VII,  we  sum¬ 
marize  our  results  while  outlining  possible  physical  realiza¬ 
tions  and  discussing  some  future  directions  of  research  in 
this  area. 


II.  PLASMON  MODES  ON  A  NANOWIRE 


The  method  for  calculating  the  electromagnetic  modes  of 
a  nanowire  is  briefly  outlined  here,  with  details  of  the  calcu¬ 
lation  given  in  Appendix  A.  We  consider  a  cylinder  of  radius 
R  of  dimensionless  electric  permittivity  e2,  which  is  centered 
along  the  z  axis  and  surrounded  by  a  second  dielectric  me¬ 
dium  ex.  While  we  are  particularly  interested  in  the  case  of  a 
conducting  nanowire  surrounded  by  some  lossless  positive 
dielectric  (Re  62<0,  ex  >0),  we  note  that  at  this  point  the 
discussion  is  quite  general.  Like  any  other  simple  geometry 
with  a  high  degree  of  symmetry,  one  can  use  separation  of 
variables  and  find  field  solutions  E,  H  to  Maxwell’s  equa¬ 
tions  in  each  dielectric  region.26,27  In  cylindrical  coordinates, 
the  electric  field  is  given  by  Ei(r)=£imEim(ki_Lp)eim(f,elk^, 
where  i=  1,2  denote  the  regions  outside  and  inside  the  cyl¬ 
inder,  respectively.  Here,  k\\  is  the  longitudinal  component  of 
the  wave  vector,  which  is  related  to  the  vacuum  wave  vector 
k0=o)/c ,  electric  permittivity  eh  and  transverse  wave  vector 
ki±  by  €jkl=k^-\-kf±,  and  m  is  an  integer  characterizing  the 
winding  of  the  mode.  The  functions  Ei  m  represent  some  nor¬ 
malized  mode  profiles.  A  similar  expression  holds  for  the 
magnetic  field  H.  For  future  reference,  we_  also  define  the 
vacuum  wavelength  \0=27r/k0,  and  kt=  \  6^0  as  the  wave 
vector  in  medium  i.  The  coefficients  £im  and  H^m  multiply¬ 
ing  the  fields  are  not  arbitrary  but,  instead,  must  satisfy  a  set 
of  equations  that  enforces  the  necessary  boundary  conditions 
at  the  dielectric  interface  p=R.  The  existence  of  a  nontrivial 
solution  requires  that  the  matrix  corresponding  to  this  linear 
system  have  zero  determinant  (detM  =  0),  which  upon  sim¬ 
plifying  yields  the  mode  equation,26,27 


m2kf\  I  1 


1  J'm(k2±R)_  1  H'm(k1±R) 

_  k2±  Jm(k2±R)  kl±Hm(k1±R)  _ 


X 


J'JhiR) 

_k21Jm(k2±R) 


k]  H'JktlR) 
k\LH,JkXLR)_  ' 


(1) 


The  above  equation,  for  example,  determines  the  allowed 
values  of  k\\  as  functions  of  m,  R,  and  er 

We  now  focus  on  the  case  of  a  subwavelength,  conducting 
metal  wire  surrounded  by  a  normal,  positive  dielectric.  In 
Fig.  2,  we  plot  the  allowed  wave  vectors  k\ |,  as  determined 
through  Eq.  (1),  as  a  function  of  R  for  a  few  lo west-order 
modes  in  m.  For  concreteness,  all  numerical  results  presented 
in  this  paper  are  for  a  silver  nanowire  (or  later,  nanotip)  at 
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FIG.  2.  (Color  online)  Allowed  plasmon  modes  k\\  as  a  function 
of  R  for  a  silver  nanowire  embedded  in  a  surrounding  dielectric 
e1=2,  for  frequency  corresponding  to  a  vacuum  wavelength  X0 
=  1  jam  and  room  temperature.  The  fundamental  (m-0)  mode,  in 
black,  exhibits  a  k\\ a  1  /R  dependence,  while  all  other  modes  are 
effectively  cut  off  as  R^O.  In  order  of  increasing  cutoff  radius,  the 
other  modes  displayed  correspond  to  \m  \  =  1 ,  2,  and  3,  respectively. 
Inset:  the  propagative  losses  for  the  fundamental  mode,  character¬ 
ized  by  the  ratio  Re£||/Imk||,  for  the  same  parameters. 

room  temperature,  \0=  1  fi m,  and  with  a  surrounding  dielec¬ 
tric  ex  =  2,  although  the  physical  processes  described  are  not 
specific  to  silver  or  to  some  narrow  frequency  range.  The 
electric  permittivity  of  the  silver  nanowire  at  this  frequency 
is  assumed  to  correspond  to  its  measured  value  in  thin  films, 
62~-50  +  0.6/.28,29  In  plotting  Fig.  2,  we  have  temporarily 
ignored  the  dissipative  imaginary  part  of  e2,  although  we  will 
address  its  effect  later.  Ignoring  Im  e2  results  in  purely  real 
values  of  kh  indicating  that  these  modes  propagate  without 
loss. 

We  first  qualitatively  discuss  the  important  features  of  the 
plasmon  modes  illustrated  in  Fig.  2  before  deriving  them 
more  carefully.  It  is  clear  from  the  figure  that  the  longitudinal 
component  of  the  wave  vector  exceeds  the  wave  vector  in 
uniform  dielectric,  h>ku  which  in  turn  causes  the  perpen¬ 
dicular  component  kil  =  ^k^-kfl=iKil  to  be  purely  imagi¬ 
nary.  Physically,  these  relationships  imply  that  the  plasmon 
modes  are  nonradiative  and  are  confined  near  the  metal/ 
dielectric  interface,  with  the  length  scale  of  transverse  con¬ 
finement  determined  by  ~  1  / Ki±.  Furthermore,  these  plas¬ 
mon  modes  cannot  couple  directly  to  radiative  fields,  which 
have  wave  vectors  k\ \  ^kx.  Of  particular  interest  is  the  behav¬ 
ior  of  the  plasmon  modes  in  the  nanowire  limit  1.  In 

this  limit,  all  higher-order  modes  (\m\  ^  1)  exhibit  a  cutoff  as 
0,  as  derived  in  Appendix  B,  while  the  m  =  0  fundamen¬ 
tal  plasmon  mode  exhibits  a  unique  k\\ &1/R  behavior.  This 
scaling  indicates  that  the  wavelength  of  these  surface  plas- 
mons  can  become  strongly  reduced  relative  to  the  free- space 
wavelength.  Physically,  in  this  limit,  the  m  =  0  mode  can  be 
interpreted  approximately  as  a  quasistatic  configuration  of 
field  and  associated  charge  density  wave  on  the  wire.  As 
such,  R  becomes  the  only  relevant  length  scale,  as  the  length 
scales  l/\ki\  associated  with  electrodynamic  behavior  be¬ 
come  unimportant.  From  the  l/R  scaling  of  k\\,  it  also  fol¬ 
lows  that  kx±^1/R.  This  implies  that  the  field  outside  the 
wire  becomes  tightly  localized  on  a  scale  <*R  around  the 
metal  surface,  leading  to  a  small  effective  transverse  mode 


area  that  scales  like  Aeff^R2.  In  particular,  it  is  possible  for 
this  mode  to  be  confined  well  below  the  diffraction  limit. 
This  effect  gives  rise  to  a  strong  interaction  strength  with 
nearby  emitters,  as  will  be  discussed  in  following  sections. 
We  note  that  this  behavior  contrasts  sharply  with  that  of,  e.g, 
a  subwavelength  normal  dielectric  waveguide  or  optical  fi¬ 
ber,  which  runs  into  a  “confinement  problem,”  where  the 
evanescent  tails  outside  the  device  become  exponentially 
large  as  /?— >  0. 30 

In  practice,  e2  is  not  purely  real  but  has  a  small  imaginary 
part  corresponding  to  losses  in  the  metal  (heating)  at  optical 
frequencies.  Its  effect  is  to  add  a  small  imaginary  component 
to  k\\,  corresponding  to  dissipation,  as  the  plasmon  propagates 
along  the  wire.  In  the  inset  of  Fig.  2,  we  plot  Re  k^/lm  k\\  for 
the  fundamental  mode  as  a  function  of  R.  This  quantity  is 
proportional  to  the  decay  length  in  units  of  the  plasmon 
wavelength  \pl  =  27r/Re  k\\.  As  R  decreases,  it  can  be  seen 
that  this  ratio  decreases  monotonically  but  approaches  a  non¬ 
zero  constant,  as  will  be  explicitly  shown  below.  For  silver  at 
\0=  1  room  temperature,  and  e1  =  2,  this  constant  is  ap¬ 
proximately  140.  The  fact  that  this  ratio  does  not  approach 
zero  even  as  R— >0  is  important  for  potential  applications,  as 
it  implies  that  the  plasmons  can  still  travel  several  times  the 
wavelength  \pl  for  devices  of  any  size.  We  also  note  that 
while  all  numbers  and  figures  presented  here  are  for  room 
temperature,  operating  at  lower  temperatures  might  some¬ 
what  reduce  the  value  of  Im  e2  due  to  decreased  losses  from 
phonon-assisted  absorption.31 

We  now  analyze  the  fundamental  mode  more  carefully. 
For  m-  0,  one  sees  in  Eq.  (1)  that  one  of  the  two  terms  on  the 
right-hand  side  must  equal  zero.  It  can  be  shown  that  setting 
the  first  term  to  zero  corresponds  to  a  TE  mode,  while  the 
other  case  corresponds  to  a  TM  mode  (see  Appendix  A).  The 
TE  mode  equation  does  not  have  any  solutions,  and  thus,  the 
fundamental  mode  is  a  TM  mode  that  satisfies  the  simplified 
equation10,32 

*2  kj  H&kl±R)_  (2) 

k2±J0(k2±R)  kl±H0(kl±R) 


The  fields  themselves  are  given  by  (see  Appendix  A) 
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ik\\k]  i  &i  | 

Ho(kl±P)P+  -jjrHo(kl±p)z 

Kx  /Cj 


jk\\z 


E?  —  b : 


“72^  J'o(k2Lp)p  +  ^jrJ0(k2±p)z 


0ik\\z 


«!  =  —ki±biH'0(kt±p)eik^4>, 

ufJ-o 


H2=—k2±b2J'0(k21p)eik^,  (3) 

(op,0 

while  the  boundary  conditions  between  the  two  dielectrics 
require  that 

bi  _  *2i  Jo(k2±R) 

b2  k1±H$cl±R)' 


035420-3 


CHANG  et  al. 


PHYSICAL  REVIEW  B  76,  035420  (2007) 


The  HR  dependence  of  k\\,  kiL  in  the  nano  wire  limit  can 
be  confirmed  mathematically  by  considering  the  nonretarded 
limit,  In  this  case,  ki±  =  \lkf-kf\  ~  ik\\  and  the  mode 

equation  (2)  reduces  to 

£2  = 

where  Im,  Km  are  modified  Bessel  functions.  The  solution  to 
Eq.  (5)  requires  k\\R  =  C  to  be  constant  and  proves  the  afore¬ 
mentioned  scaling  law  for  k\\.  It  is  also  straightforward  to  see 
that  when  e2  acquires  a  small  imaginary  component,  the  con¬ 
stant  C  becomes  complex  as  well,  and  that  Refc||/Imfc|| 
=Re  C/Im  C  takes  on  some  fixed,  nonzero  value.  No  closed- 
form  solution  exists  for  the  equation  above,  although  when 
|&I|R|<^1  (corresponding  to  large  \e2l £\ |),  the  equation  as¬ 
ymptotically  approaches 


ei  (y-  log  2  +  log  C)(C)2  ’ 

where  y~  0.577  is  Euler’s  constant. 

Finally,  it  should  be  noted  that  the  components  of  E*  in 
Eq.  (3)  are  proportional  to  k\\ki±  or  kf±,  while  H*  is  propor¬ 
tional  to  ki±.  Thus,  in  the  nanowire  limit  when  &u,|A7-±| 
ocl  /R,  the  magnetic  fields  are  a  factor  of  R  smaller  than  the 
electric  fields,  which  is  consistent  with  this  mode  being 
roughly  a  quasistatic  configuration. 

III.  SPONTANEOUS  EMISSION 
NEAR  A  METAL  NANO  WIRE 

The  small  mode  volume  associated  with  the  fundamental 
plasmon  mode  of  a  nano  wire  offers  a  possible  mechanism  to 
achieve  strong  coupling  with  nearby  optical  emitters,  in  anal¬ 
ogy  to  the  methods  of  Refs.  22-25.  In  this  section,  we  derive 
more  rigorously  the  interaction  between  an  emitter  and  an 
infinite,  cylindrical  nanowire,  and  show  that  under  certain 
circumstances  the  small  mode  volume  leads  to  strongly  pref¬ 
erential  spontaneous  emission  into  the  guided  plasmon 
modes  via  a  mechanism  equivalent  to  the  Purcell  effect33  in 
cavity  QED. 

The  spontaneous  emission  rate  of  a  dipole  emitter,  in  gen¬ 
eral,  becomes  altered  from  its  free-space  value  in  the  pres¬ 
ence  of  some  dielectric  body.  In  our  system  of  interest,  the 
dipole  can  possibly  lose  power  radiatively  to  propagating 
photon  modes  through  excitation  of  the  guided  plasmon 
modes  or  through  nonradiative  loss  (heating)  in  the  wire.  The 
dipole  in  consideration  can  physically  be  formed  by  a  single 
atom,  a  defect  in  a  solid-state  system,  or  any  other  system 
with  a  dipole-allowed  transition.  In  Sec.  Ill  A,  we  calculate 
the  radiative  and  nonradiative  rates  using  a  quasistatic  ap¬ 
proach,  while  waiting  until  Sec.  Ill  B  to  treat  the  plasmon 
decay  rate  more  thoroughly.  In  Sec.  Ill  C,  we  show  how  the 
efficiency  of  emission  into  the  plasmon  modes  can  be  opti¬ 
mized  to  yield  Purcell  factors  in  excess  of  ~103,  and  discuss 
the  physical  origins  of  this  limit. 

A.  Radiative  and  nonradiative  decay  rates 

In  this  section,  we  derive  formulas  for  the  decay  rates  of  a 
dipole  near  a  metal  nanowire  into  radiative  and  nonradiative 


d  =(v'2-Vq)/ 2 


FIG.  3.  (Color  online)  (a)  A  dipole  emitter  positioned  a  distance 
d  from  the  center  of  a  nanowire.  For  the  calculations  presented,  the 
orientation  of  the  dipole  moment  is  assumed  to  be  along  p.  (b)  A 
dipole  emitter  positioned  near  the  end  of  a  nanotip.  In  terms  of 
parabolic  coordinates,  the  physical  distance  d  between  the  end  of 
the  tip  and  the  emitter  is  given  by  d={v,2-v2^)l2,  while  the  curva¬ 
ture  parameter  w  =  2v$.  For  the  nanotip,  only  a  dipole  oriented  par¬ 
allel  to  z  couples  to  the  fundamental  plasmon  mode. 

channels.  This  calculation  closely  follows  that  of  Ref.  34,  but 
is  briefly  presented  here  for  completeness. 

It  is  well  known  that  spontaneous  emission  rates  can  be 
obtained  via  classical  calculations  of  the  fields  due  to  an 
oscillating  dipole  near  the  dielectric  body,35  and  this  method 
will  be  employed  here.  Specifically,  we  consider  an  (classi¬ 
cal)  oscillating  dipole  p 0e~la)t  positioned  a  distance  d  from 
the  center  of  the  wire,  and  wish  to  calculate  the  total  fields  of 
the  system  [see  Fig.  3(a)].  Before  proceeding  further,  we  first 
note  that,  in  general,  the  solution  will  depend  on  the  orien¬ 
tation  of  the  dipole  moment  in  addition  to  its  position.  Again 
taking  £  to  be  along  the  axis  of  the  wire,  we  note  that  the 
fundamental  surface  plasmon  mode  will  not  couple  to  a  di¬ 
pole  moment  p0^  0  oriented  along  the  azimuthal  axis,  due  to 
the  azimuthal  symmetry  of  this  mode.  In  the  following  cal¬ 
culations  for  the  nanowire,  for  concreteness  we  will  consider 
a  dipole  oriented  along  the  radial  direction  (po  ap)  while 
noting  that  the  calculation  for  an  orientation  along  £  yields 
qualitatively  similar  results.  For  nanostructures,  one  can 
make  a  simplification  by  considering  the  fields  in  the  quasi- 


static  limit  (H  =  0),34  which  satisfy 

V  D  =  pext. 

(7) 

V  X  E  =  0. 

(8) 

Here,  Pext(r)  is  the  external  charge  density.  In  taking  the 
quasistatic  limit,  we  essentially  assume  that  we  are  consider¬ 
ing  length  scales  r  over  which  the  phases  ktr  associated  with 
the  electrodynamic  Green’s  function  G*(r)  ~  elk'r / r  are  neg¬ 
ligible,  so  that  this  function  can  be  approximated  by  the  elec¬ 
trostatic  version,  G;(r)~l/r.  In  the  system  of  interest,  the 
external  source  is  a  dipole  located  at  position  r'  outside  the 
wire  (with  radial  coordinate  p'=d),  so  that 

Pext(r,r')  =  (Po-V')^(r-r').  (9) 

For  simplicity,  we  omit  the  harmonic  time  dependence  from 
our  expressions  for  the  source  and  all  fields.  Note  that  the 
£(r-r')  term  above  corresponds  to  a  (unitless)  point  charge 
source,  while  the  operator  (p0-V')  generally  converts  the 
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point  charge  solution  to  that  of  a  dipole.  It  is  therefore  con¬ 
venient  to  write  E*  in  a  similar  form, 


Ej(r,r')  =  -  V(p0  •  V')<P,(r,r'),  (10) 

where  <3Y(r,r')  are  “pseudopotentials”  that  satisfy  V2^ 
=-S( r-r')/e0e1  and  V2<E2  =  0.  Here,  the  indices  1  and  2 
again  denote  the  regions  outside  and  inside  the  cylinder,  re¬ 
spectively.  physically  corresponds  to  the  potential  due  to  a 
point  charge  at  r',  while  the  dipole  potential  follows  from 
®<«p=(Po-V')®. 

To  solve  for  the  fields,  it  is  convenient  to  further  separate 
dfi  into  “free”  and  “reflected”  components  O0  and  Or,  re¬ 
spectively,  where  Or  represents  a  source-free  contribution 
that  ensures  that  boundary  conditions  are  satisfied,  and  O0  is 
the  solution  for  a  point  charge  in  a  medium  of  uniform  elec¬ 
tric  permittivity  ex .  We  will  expand  the  known  source  term 
<E0  in  a  basis  appropriate  for  the  cylindrical  geometry  and 
expand  the  source-free  terms  Or  .  2  in  a  similar  basis  that  sat¬ 
isfies  Laplace’s  equation  (V2Or2  =  0).  The  unknown  coeffi¬ 
cients  multiplying  the  basis  functions  of  2  will  then  be 
determined  by  enforcing  the  proper  boundary  conditions  at 
the  dielectric  interface.  These  expansions  are  given  by 


%(ry)  = 


_L _ l_ 

47T€0el  |r  —  r ' | 

00 

:  — -j - 2  (2  -  Sm0)cos[m(</>-  </>')] 

Z7T  6061  m=0 
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X  I  dh  cos[/z(z  -  z')]Km(hp')Im{hp ) 
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X  |  dham(h)cos[h(z- z')]Km(hp')Km(hp), 

(12) 


a>2(r’r')  =  7^X  ^2~  <WcosM0~  ^,)1 


"0  m=0 


roo 

l 


X  I  dh/3m(h)cos[h(z-z')]Km(hp')Im(hp), 

(13) 


where  am(/z)  and  /3m(h)  thus  far  are  unknown  amplitude  co¬ 
efficients.  We  obtain  a  set  of  two  coupled  equations  for 
am{h)  and  by  requiring  continuity  of  and  Dx  at  the 

boundary,  p=R.  Because  of  the  translational  symmetry  of  the 
system,  these  equations  are  uncoupled  in  h  and  can  easily  be 
solved  (this  is  in  contrast  to  the  case  where  translational 
symmetry  is  broken  due  to  surface  roughness,  as  discussed  in 
Sec.  VI).  The  solutions  are  given  by34 


^m(^) 


(6-  1  )l'm(hR)Im(hR) 
eJmm)K'm{hR)  -  e2Km(hR)rm(hR) 9 


UhR)K'JhR)-Km(hR)rm(hR) 
€lIm(hR)K'm(hR)  -  €2Km(hR)I'm(hR)  ’ 


(14) 


where  we  have  defined  e=e2lex.  Note  that  Eq.  (14)  along 
with  Eqs.  (1 1)— (13)  gives  the  total  electric  field  of  the  dipole 
and  wire  system. 

To  calculate  the  radiative  emission  into  free  space,  we 
consider  the  far-field  properties  of  the  system.  Physically,  the 
presence  of  the  emitter  induces  some  dipole  moment  £p  in 
the  nanowire,  which  results  in  a  total  radiated  power  propor¬ 
tional  to  the  square  of  the  total  dipole  moment  of  the  system, 
|p0+ ^p|2.  We  can  determine  <5p  by  finding  the 
dipolelike  contribution  to  the  reflected  potential  d>^r 
=  (po *  VA)Or(r ,  r^)  far  away  from  the  source,  which  on  physi¬ 
cal  grounds  must  behave  like  p-2  for  large  p.  It  is  straight¬ 
forward  to  show  that  the  m=  1  term  in  Eq.  (12)  is  responsible 
for  this  contribution,  with  all  other  m  terms  yielding  faster 
decays  in  p.  Because  of  the  asymptotic  behavior  of  Km(x) 
~  V 7r/2xe~x  when  1,  for  large  p  the  integrand  in  Eq.  (12) 

is  significant  only  over  a  small  region  h  ^  p-1.  As  a  result,  we 
can  safely  replace  Kx(hpr)  and  ax(h)  by  their  expansions 
around  h-  0.  After  this  simplification,  the  integral  can,  in 
fact,  be  evaluated  exactly  and  yields 


<f> 


(m=l)  . 


4776^6+1  C0S(<^'V  [p2  +  (Z-Z’ )2]3/2’ 

(15) 


with  a  corresponding  reflected  potential 
^r1)(r,r')  =  (Po-V')^m=1)(r,r') 

V)r\  £  —  1 

■  COs(<^>-  4>  )  2  ,  2  ,  „2\3/2 


Po  1  . R2 


4T7e0€i  6+1 


d2(p2  +  z2)m ’ 

(16) 


for  the  choice  of  parameters  p'  =d,  z'=  0,  and  p0=p0p.  Com¬ 
paring  Eq.  (16)  to  the  potential  due  to  a  dipole  <5p  in  uniform 
dielectric  eb  ^=4^>'Jp,  we  can  readily  identify 


<5p=p0 


6-1  j?2 
6+  1  d2 


(17) 


as  the  induced  dipole  moment  in  the  wire,  from  which  it 
follows  that  the  radiative  spontaneous  emission  rate  is  given 

by34 


1  + 


6-  l/?2 
6+  1  d2 


2 


(d^R). 


(18) 


Here,  ro  is  defined  to  be  the  spontaneous  emission  rate  of  the 
emitter  in  uniform  dielectric  61.36  Away  from  the  plasmon 
resonance  (6~-l),  the  radiative  decay  rate  changes  slightly 
from  T o  and  reflects  some  moderate  change  in  the  radiative 
density  of  states  in  the  vicinity  of  the  nano  wire. 

To  calculate  the  other  decay  rates,  one  utilizes  the  fact  that 
the  total  power  loss  of  an  oscillating  dipole  is  proportional  to 
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the  electric  field  in  quadrature  at  the  dipole’s  location,  spe¬ 
cifically,  rto^rclmfpo-E^r'  ,r')].  Having  divided  up  E{  into 
free  and  reflected  components,  the  contribution  to  Ex  from 
the  free  field  simply  is  associated  with  the  decay  rate  in 
uniform  dielectric  and  thus,  we  concentrate  on  the  con¬ 
tribution  from  Or(r,r').  First,  we  note  that  the  coefficient 
a0(h)  derived  in  Eq.  (14)  contains  a  pole  at  the  point  where 
the  denominator  vanishes.  This  pole  corresponds  to  an  exci¬ 
tation  of  a  natural  mode  (the  fundamental  plasmon  mode)  of 
the  system.  This  can  immediately  be  seen  by  comparing  the 
denominator  of  a0  to  Eq.  (5),  which  gives  the  plasmon  mode 
in  the  nano  wire  limit.  The  pole  lies  at  h  =  C/R  and  agrees 
with  the  plasmon  wave  vector  derived  in  Sec.  II,  as  expected. 
Evaluating  the  contribution  of  this  pole  to  Er(r',r')  gives 
the  decay  rate  into  the  fundamental  plasmon  mode,  and  is 
discussed  more  carefully  in  the  next  section.  At  the  same 
time,  in  the  limit  d^R,  one  expects  some  type  of  divergence 
to  occur  in  the  nonradiative  decay  rate.  Physically,  as  the 
emitter  approaches  the  wire  edge,  the  divergence  results 
from  the  large  currents  in  the  wire  generated  by  the  near  field 
of  the  dipole  and  their  resulting  dissipation.  We  can  find  the 
leading-order  term  to  this  divergent  decay  rate  by  carefully 
evaluating  the  leading-order  divergence  in  the  reflected  field. 

The  mathematical  origin  of  the  divergence  as  R  is  the 
significant  contribution  to  <Fr  of  an  infinite  number  of  terms 
in  m.  Specifically,  in  this  limit,  for  a  dipole  oriented  along  p, 


■  non-rad 
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dhfm(h,d,R).  (19) 


The  asymptotic  behavior  of  the  functions  fm  is  given  by 


fm(h,d,R) 
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(20) 


From  the  above  expressions,  we  see  that  fm  as  a  function  of 
h  has  a  characteristic  width  of  about  [2 (d-R)]-1,  yet,  at  the 
same  time,  the  quantity  m(R/d)2m  reaches  a  maximum 
around  m  ~  2{R-d)  as  This  confirms  the  nonvanishing 

contribution  of  an  infinite  number  of  terms  in  m  to  the  decay 
rate.  The  exact  behavior  of  the  functions  fm  is  well  modeled 
by  a  Lorentzian  approximation, 


fm(h,d,R) 
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ld2ex 
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(21) 


which  allows  the  integration  and  sum  in  Eq.  (19)  to  be  per¬ 
formed  exactly.  The  resulting  decay  rate  is  given  by 


r 


nonrad 


r0 


16 kl(d-R)34n  Im(  e  +  1 ) 


(22) 


Note  that  for  |e|>l  and  small  Im  6,  Im(^j-) 

~2  Im  6/ (Re  e)2,  which  makes  it  clear  that  the  nonradiative 
spontaneous  emission  rate  is  proportional  to  the  dissipative 
part  of  the  electric  permittivity. 


B.  Decay  rate  into  plasmon  modes 


In  this  section,  we  quantify  the  spontaneous  emission  rate 
Tpl  of  a  dipole  into  the  surface  plasmon  modes  on  a  nanowire 
by  evaluating  the  contribution  of  the  pole  in  a0(h)  to  the  field 
seen  by  the  dipole.  Before  proceeding  further,  we  first  note 
that  in  the  presence  of  metal  losses,  the  distinction  between 
Tpl  and  Tnonrad  is  not  perfectly  well  defined,  since  the  plas- 
mons  eventually  dissipate  due  to  heating  as  well.  Thus,  for 
concreteness,  we  will  define  Tpl  to  be  the  decay  rate  resulting 
from  the  pole  in  the  limit  that  Im  62  =  0,  and  take  the  plasmon 
wave  vector  k\\  and  C  to  be  purely  real  in  this  section.  In 
particular,  for  a  dipole  oriented  along  p, 
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,  (23) 


pole 


where  we  have  explicitly  indicated  that  we  are  interested  in 
the  pole  contribution  to  the  expressions  above.  The  behavior 
of  a0(h)  around  the  pole’s  vicinity  is  given  by 


a0(h)  = 


1  (ei-eMO  /0(O 


(h  -  C/R)R 


dx(C)  ’ 
dx 


(24) 


where 


X(x)  =  e\I{)(x)K’(i(x)  -  e2K0(x)I'0(x) .  (25) 

Substituting  this  into  Eq.  (23)  then  yields  the  decay  rate 


^  K](Cd!R) 
lpl=apll°  (k0R)3 


K\(Klld) 
"P'r°  (k0R)3 


(26) 


where  we  have  identified  kx±^C/R  in  the  nano  wire  limit. 
The  coefficient  apl  is  given  by 


«pi  = 


3(ft  ~  g2)  C2/1(C)/0(C) 

6[/2  dx(C)/dx 


(27) 


and  most  importantly  depends  only  on  el  2. 
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The  ~(1  /R)3  scaling  of  the  spontaneous  emission  rate 
into  plasmons  indicates  that  this  process  can  be  strongly  en¬ 
hanced  by  using  a  small  wire.  One  can  gain  some  physical 
understanding  of  the  result  and  its  relation  to  the  Purcell 
effect  in  cavity  QED  by  considering  the  quantum  mechanical 
derivation  of  the  spontaneous  emission  rate  via  Fermi’s 
golden  rule.  This  rule  states  that,  once  the  plasmon  modes 
are  quantized,  the  decay  rate  is  given  by 

rpi  =  2i7g2(r,  «)£>(«),  (28) 

where  g(r,co)  is  the  position-dependent,  dipole-field  interac¬ 
tion  matrix  element,  and  D(a))  is  the  plasmon  density  of 

states  on  the  nano  wire.  Here,  g(r ,  co) = $?  -  £(r)  /  h  itself  de¬ 
pends  on  the  dipole  matrix  element  p  =  (e\er\g)  of  the  emit¬ 
ter,  where  | g)  (\e))  is  the  ground  (excited)  state  involved  in 

the  transition,  and  on  the  electric  field  per  photon  £(r)  for  the 
surface  plasmon  modes.  On  dimensional  grounds,  the  elec- 
tric  field  per  photon  can  be  written  in  the  form  £(r) 
=  ^ha)/ e0Veff[E(r)/Emax\,  where  Veff  is  an  effective  mode 
volume  characterizing  the  confinement  of  each  mode,  E(r)  is 
the  classical  field  profile  for  the  surface  plasmon  modes 
given  in  Eq.  (3),  and  £WflJC=max|E(r)|.  The  mode  volume  can 
be  estimated  by  normalizing  the  field  energy  to  one  quantum 
(again,  ignoring  Im  62), 

h(o~L  j  d2pe0^[«e(p,«)]|£(p)|2.  (29) 


Here,  the  two-dimensional  integral  is  performed  over  the  di¬ 
rections  transverse  to  the  wire  axis,  and  L  is  the  quantization 
length,  set  by  the  wire  length,  which  is  assumed  to  be  much 
longer  than  all  other  relevant  length  scales.  In  the  end,  L  will 
disappear  from  the  physical  quantities  of  interest  (e.g.,  the 
spontaneous  emission  rate).  The  integrand  appearing  in  Eq. 
(29)  gives  the  correct  expression  for  the  classical  electric 
energy  density  for  a  field  in  a  dispersive  medium.27  Because 
the  fields  are  primarily  electric  in  the  nanowire  limit,  we 
have  safely  ignored  the  magnetic  field  contribution  to  the 
energy  density.  Performing  this  normalization  for  the  funda¬ 
mental  plasmon  modes,  one  finds  that  Veff  <*R2L,  which 
physically  implies  that  these  modes  are  transversely  confined 
to  an  area  on  the  order  of  the  wire  size.  The  constant  of 
proportionality  depends  only  on  the  electric  permittivities  of 
the  system.  For  a  dipole  oriented  along  p,  one  then  finds  that 
g2(r)^Kl(Klld)/(R2L),  which  shows  that  the  strong  field- 
emitter  coupling  directly  arises  from  the  tight  field  confine¬ 
ment.  At  the  same  time,  the  density  of  states  for  the  plasmon 
modes  is  given  by  D(a))  ~  (L/2Tr)(dk\\/doj),  and  is  inversely 
proportional  to  the  group  velocity.  The  group  velocity  for 
plasmons  is  strongly  reduced  due  to  the  large  plasmon  wave 
vector  k\\  1  /R  on  a  nanowire,  and  consequently,  this  gives 

rise  to  a  — HR  enhancement  in  the  density  of  states.  Com¬ 
bining  all  of  these  results  into  Eq.  (28)  yields 


rPi 


rn 


^iUhqQ2 

(M)3 


(30) 


where  again  the  constant  of  proportionality  depends  only  on 
the  electric  permittivities.  We  emphasize  that  in  the  quantum 


mechanical  interpretation,  one  achieves  a  factor  of  ~1  /R 
enhancement  in  the  spontaneous  emission  rate  due  to  a  re¬ 
duction  in  the  plasmon  group  velocity  and  a  ~  1  /R2  enhance¬ 
ment  due  to  the  small  effective  mode  area. 


C.  Purcell  factor  of  a  nanowire 

Comparing  the  spontaneous  emission  rates  given  by  Eqs. 
(18),  (22),  and  (26),  we  now  qualitatively  discuss  the  behav¬ 
ior  one  should  expect  as  the  position  of  the  emitter  is  varied. 
In  the  limit  that  d/R>  1,  the  emitter  feels  no  effect  from  the 
wire  and  the  total  spontaneous  emission  rate  is  close  to  the 
radiative  rate  ro  in  a  uniform  dielectric  medium.  As  one 
brings  the  emitter  closer  to  the  wire  surface,  the  change  in 
the  electromagnetic  mode  structure  near  the  wire  results  in 
some  modified  radiative  decay  rate  rrad  which  never  exceeds 
approximately  4ro  for  large  |e|.  When  the  emitter  position  d 
approaches  d~l/  ki±~R/\C\9  the  emitter  starts  to  interact 
with  the  localized  plasmon  fields,  with  a  corresponding  rate 
of  emission  into  plasmons  scaling  with  wire  size  like  l/R3. 
The  emission  rate  into  plasmons  continues  to  grow  as  the 
emitter  is  brought  even  closer  to  the  wire  edge,  d^R.  How¬ 
ever,  the  efficiency  or  probability  of  plasmon  excitation 
eventually  decreases  due  to  the  large  nonradiative  decay  rate 
experienced  by  the  dipole  very  near  the  wire,  which  diverges 
like  1  /(d-R)3.  We  thus  expect  some  optimal  efficiency  of 
spontaneous  emission  into  the  plasmon  modes  to  occur  when 
the  emitter  is  positioned  at  a  distance  0(R)  away  from  the 
wire  edge,  and  for  this  optimal  efficiency  to  improve  as  R 
->0. 

The  efficiency  of  coupling  to  the  plasmon  modes  can  be 
characterized  by  a  “Purcell  factor”  that  is  defined  by  the  ratio 
P=rpl/r',  where  r'  =  rrad+rnonrad  denotes  the  total  emission 
rate  into  channels  other  than  the  fundamental  plasmon  mode. 
Generally,  P  depends  on  both  the  wire  size  and  the  position 
d  of  the  emitter.  In  Fig.  4(a),  we  have  numerically  evaluated 
the  spontaneous  emission  rates  and  plotted  the  Purcell  factor 
as  a  function  of  wire  radius  R ,  optimized  over  the  emitter 
position.  As  R—> 0,  the  optimized  Purcell  factor  P  exceeds 
~103,  indicating  that  the  probability  of  emission  into  the 
plasmons  approaches  almost  unity.37  Examining  this  limit 
more  carefully,  it  can  be  shown  that  the  “error  rate,”  in  fact, 
approaches  a  small  quantity  r'/rtotaloclm  6/ (Re  e)2,  which 
explicitly  indicates  that  this  process  is  ultimately  limited  by 
material  losses.  Again,  we  emphasize  that  these  properties 
are  specifically  a  result  of  the  conducting  properties  of  the 
nano  wire.  This  can  be  contrasted  with  emission  into  the 
guided  modes  of  a  subwavelength  optical  fiber,  which  drops 
exponentially  as  R  — >  0  due  to  the  weak  confinement  of  these 
guided  modes.34 

Finally,  we  note  that  while  the  decay  rates  obtained  in  the 
quasistatic  approximation  provide  a  simple  understanding  of 
the  system,  the  quantitative  predictions  may  begin  to  break 
down  even  for  moderate- sized  wires  (when  \kj\R^l).  Physi¬ 
cally,  the  effective  mode  area  Aeff  and  group  velocity  do)/dk\\ 
start  to  scale  more  slowly  than  ~R2  and  ~R,  respectively, 
for  moderate- sized  wires.  The  quasistatic  approximation  can 
then  significantly  underestimate  the  decay  rate  into  plas¬ 
mons,  rpl,  and  the  corresponding  Purcell  factors.  This  may 
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FIG.  4.  (Color  online)  (a)  Solid  line:  Maximum  Purcell  factor  P=rpl/r'  for  a  nano  wire,  plotted  as  a  function  of  R  and  optimized  over 
the  emitter  position.  Here,  P  is  calculated  in  the  quasistatic  approximation.  Dashed  line:  Same  quantity,  obtained  by  exact  electrodynamic 
calculations.  Dotted  line:  Effective  Purcell  factor  P(R)  =  TvX(R)  IT(R)  for  a  nanotip  of  final  radius  R ,  optimized  over  nanotip  shape  and  emitter 
position.  Solid  points:  Same  quantity,  calculated  using  numerical  simulations  (boundary  element  method),  (b)  Contour  plot  of  log10  P  for  a 
nanowire,  as  functions  of  R  and  d/R.  (c)  Contour  plot  of  log10  P  for  a  nanotip,  as  functions  of  w  and  dlw. 


be  important  experimentally,  for  example,  since  it  implies 
that  one  need  not  resort  to  very  small  wires  to  achieve  mod¬ 
erate  Purcell  factors.  Thus,  using  the  electrodynamic  Green’s 
function  methods  of  Refs.  34  and  38,  we  have  also  obtained 
exact  electrodynamic  solutions  of  the  fields  created  by  an 
oscillating  dipole  near  a  cylinder,  and  used  these  solutions  to 
find  the  corresponding  decay  rates.  In  Fig.  4(a),  we  have 
plotted  the  numerically  optimized  Purcell  factors  using  these 
exact  electrodynamic  results.  It  can  be  seen  that  the  quasi¬ 
static  and  exact  results  agree  well  in  the  nanowire  limit,  but 
for  larger  wires,  the  exact  results  predict  significantly  larger 
values  of  P.  In  Fig.  4(b),  we  plot  the  Purcell  factor  as  func¬ 
tions  of  both  wire  size  and  emitter  position  (using  the  exact 
results).  Here,  it  can  be  seen  that  achieving  a  large  Purcell 
factor  does  not  depend  too  sensitively  on  the  placement  of 
the  emitter,  which  indicates  that  the  strong  coupling  is  a  ro¬ 
bust  effect. 

IV.  SPONTANEOUS  EMISSION  NEAR  A  NANOTIP 

In  Secs.  II  and  III,  we  derived  and  discussed  the  physics 
of  plasmon  modes  on  a  nanowire  and  the  spontaneous  emis¬ 
sion  of  a  nearby  dipole  emitter.  It  was  seen  that  the  tight 
transverse  confinement  of  the  plasmon  modes  leads  to  a  large 
Purcell  factor  for  an  optimally  positioned  dipole  emitter  as 
the  wire  size  decreases.  At  the  same  time,  however,  it  is 
evident  that  the  tighter  confinement  is  accompanied  by  en¬ 
hanced  losses  as  the  plasmon  propagates,  and  also  by  a  large 
reduction  in  the  plasmon  wavelength  Xpl  that  could  make 
outcoupling  more  difficult.  Such  factors  could  clearly  impose 
limits  for  applications  such  as  quantum  information,  but  can 
be  circumvented  with  simple  design  improvements.  In  this 
section,  we  investigate  one  specific  design,  a  metallic  nan¬ 
otip  [see  Fig.  3(b)].  As  in  the  nano  wire  case,  one  expects  a 
subwavelength  plasmon  mode  volume,  determined  in  this 
case  by  the  tip  curvature,  and  an  associated  enhancement  of 
emission  into  the  plasmon  modes.  Here,  however,  one  ex¬ 
pects  propagative  losses  to  be  less  severe,  because  the  plas- 
mons  will  become  less  tightly  confined  as  they  propagate 
along  the  tip  and  as  the  tip  expands.  This  expansion  also 
results  in  a  gradual  increase  in  the  plasmon  wavelength  along 
the  nanotip.  For  the  nanotip,  we  are  not  able  to  obtain  full 


electrodynamic  solutions  for  the  plasmon  modes.  However, 
in  a  manner  similar  to  that  described  in  Sec.  Ill,  we  will 
calculate  all  of  the  relevant  decay  rates  in  the  quasistatic 
limit  and  describe  an  approximate  method  to  calculate  the 
propagative  losses  along  the  nanotip.  We  will  also  compare 
these  results  to  those  obtained  via  fully  electrodynamic  nu¬ 
merical  simulations,  and  we  find  that  these  two  approaches 
agree  closely. 

In  the  following,  we  will  consider  a  nanotip  whose  sur¬ 
face  can  be  parametrized  as  a  paraboloid  of  revolution  with 
symmetry  along  the  z  axis  [see  Fig.  3(b)].  Specifically,  we 
suppose  that  the  surface  of  the  nanotip  is  described  by 


z 


1 

2 


2  2 

x+y 


(31) 


a  paraboloid  of  revolution  with  apex  at  z--v\l2  (the  reason 
for  the  offset  of  the  apex  will  become  apparent  below).  We 
now  introduce  a  transformation  to  parabolic  coordinates, 


x-uv  cos  </>, 

(32) 

y  =  uv  sin  </>, 

(33) 

z=^(u2-v2). 

(34) 

While  these  coordinates  may  seem  awkward  (note,  for  ex¬ 
ample,  that  u  and  v  have  units  of  \length),  they  are  conve¬ 
nient  for  deriving  expressions  for  the  fields  and  spontaneous 
emission  rates,  which  we  will  then  express  in  more  “natural” 
coordinates  at  the  end  of  the  calculation.  In  parabolic  coor¬ 
dinates,  the  nanotip  profile  of  Eq.  (31)  is  defined  by  a  surface 
of  constant  v  =  v0.  More  generally,  any  constant  v  defines 
some  paraboloid  of  revolution  in  this  system,  while  the  unit 
vectors  u  and  v  run  normally  and  tangentially  to  these  sur¬ 
faces,  respectively. 

Now,  as  in  the  nanowire  case,  we  are  interested  in  finding 
the  quasistatic  field  solution  for  a  point  charge  source  in  the 
vicinity  of  the  nanotip,  from  which  we  can  obtain  the  field 
due  to  a  dipole  p0.  In  particular,  we  seek  solutions  of  the 
total  field  of  the  form  Eq.  (10)  with  appropriate  boundary 
conditions.  Like  before,  we  separate  the  pseudopotential 
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outside  the  nanotip  into  its  free  and  reflected  components 
O0r,  and  use  an  integral  representation  of  the  free  pseudo¬ 
potential  suitable  for  parabolic  coordinates, 

1  f00 

<I>o(r,r')  =  - - 2  (2-  Smfi)cosm((f>-  (f>')  dq 

Z7Te0elm=0  Jo 

XqJm(qu)Jm(qu’ )Im(qv)Kjqv’ )  (v  <  v'). 

(35) 


2  r  oo  3 

Po  •  Er(r',r')  =  -  —  dq^a0{q)K\{qv')  (p«  II  £), 

277e0eiJ0  v 

(40) 

as  shown  in  Appendix  C.  Because  of  the  similarity  of  the 
decay  rate  calculations  with  those  in  Sec.  Ill,  we  simply  state 
the  results  here,  while  providing  more  details  in  Appendix  C. 
For  a  dipole  positioned  along  the  z  axis  at  v  =  v ' ,  the  radiative 
and  nonradiative  spontaneous  emission  rates  are  given  by 


Here,  the  primed  coordinates  denote  the  position  of  the  di¬ 
pole.  Because  O0  fully  accounts  for  the  point  source,  Or2 
then  satisfy  Laplace’s  equation.  Using  separation  of  vari¬ 
ables,  it  is  straightforward  to  show  that  the  solutions  to 
Laplace’s  equation  are  given  in  parabolic  coordinates 
by  ~Jm(qLl)Gl,Jqv)eim4>,  where  Ghm(qv)=Km(qv)  and 
G2 ,m(clv)=Im(clv)  are  nondivergent  functions  in  their  regions 
of  applicability.  We  then  define  the  following  expansions: 

1  f00 

^r(r>r')  =  Z - 2  (2-  Smfi)cosm(<f>-  <f>')  dq 

27r606l  m=0  Jo 

Xqam(q)Jm(qu’  )Kjqv’  )Jm(qu)Km(qv) , 

(36) 


Efad  _  .  ^1  ^2  ^0 

T-I  —  1  F  f  9 

Tq  +  62  V 


(Pollz), 

(Po  -L  z). 


(41) 


and 


r 


nonrad 


r0 


3 _ 1 

16klel/2  v'3(v'  -  l>o)3 


r  3  i 

1  nonrad  __  _ _ _ _ 

r0  &k30el/2  v,3(v'  -v0)3 


flZll) 

62+  61/ 


(Po  -L  z), 

(Pollz)- 

(42) 


1  f  °° 

<l)2(r,r')  =  7 - 2  (2  -  4n,o)cos  m(4>-  4>')  dq 

l7Te0e\  m=0  Jo 

Xq/3m(q)Jm(qu')Km(qv')Jm(qu)IJqv),  (37) 

where  the  coefficients  a  and  /J  are  determined  by  imposing 
boundary  conditions  at  the  nanotip  surface  v  =  v0.  Enforcing 
the  continuity  of  O  and  D±  yields 


am(q)  = 


(ji  ~  g2)4(g^o)4(^o) 

^Xn(qvo)Km(qv0)  -  e1Im(qv0)K'm(qv0)  ’ 


(38) 


Pm(q)  = 


-  ei[Im(qvo)K'm(qv0)  -  l'm(qvo)Km(qv0)] 
eil'miqvoXJqvo)  -  eiIm(qv0)K'm(qv0)  ' 


Note  that  the  coefficients  am(q ),  along  with  Eq.  (36),  com¬ 
pletely  determine  the  reflected  field. 

The  calculation  of  the  radiative  and  nonradiative  sponta¬ 
neous  emission  rates  proceeds  in  the  same  manner  as  the 
nanowire  case.  To  calculate  rrad,  we  again  look  for  a  dipole 
term  in  the  far  field  (large  v)  that  corresponds  to  an  induced 
dipole  moment  £p  in  the  nanotip,  and  then  use  the  relation¬ 
ship  rradoc|p0+<Sp|2.  At  the  same  time,  we  look  for  a  diver¬ 
gent  contribution  to  the  reflected  field  at  the  dipole  location 
as  its  position  v'  approaches  v0,  which  yields  the  leading 
term  of  the  nonradiative  decay  rate  through  rnonrad 
ocIm[p0-Er(r'  ,r')].  For  concreteness,  we  will  consider  a  di¬ 
pole  that  is  located  along  the  z  axis  [V=0,  also  see  Fig. 
3(b)],  which  yields 


Po-Er(r',r')  =  - 


2 

Po  I 

477e06i  Jo 


dq—2ax{q)K\(qv')  (p0±£). 


Finally,  we  consider  the  decay  rate  into  the  fundamental 
plasmon  mode  of  the  nanotip,  which  is  associated  with  the 
contribution  of  the  poles  in  the  integrand  of  Eq.  (40)  to 
Im[p0-Er(r'  ,r')].  Examining  the  solutions  to  a0  l  given  in 
Eq.  (38),  one  finds  that  a1  has  no  pole  in  the  range 
s^oo.  Physically,  the  absence  of  a  pole  means  that  a  dipole 
simultaneously  oriented  perpendicular  to  z  and  located  along 
the  z  axis  does  not  excite  the  fundamental  plasmon  mode  of 
the  nanotip.  This  is  easily  understood  since  a  dipole  oriented 
this  way  is  antisymmetric  with  respect  to  180°  rotations 
about  £,  while  the  plasmon  mode  is  symmetric.  On  the  other 
hand,  a0  does  have  a  pole  corresponding  to  plasmon  excita¬ 
tion.  This  pole  is  located  at  q0=C/v0 ,  where  C  is  the  solution 
to  Eq.  (5).  q0  can  be  considered  to  be  the  “wave  vector”  of 
the  fundamental  mode  (albeit  in  parabolic  coordinates). 
Evaluating  the  contribution  of  this  pole  to  the  field  is 
straightforward  and  yields  a  plasmon  decay  rate 


rpl  377  c3  2  .M-ejimioiQ 
r 0~k34,2vy2  l{q°v  ’  dx(c)/dx 


(Polls). 

(43) 


7^  =  0  (Po  -L  z) ,  (44) 

1  0 

where  x(x)  is  defined  in  Eq.  (25).  As  in  the  nano  wire  case, 
the  decay  rate  rpl  into  the  plasmon  mode  given  by  Eq.  (43)  is 
evaluated  in  the  limit  that  Im  62  =  0,  such  that  q0  and  C  are 
purely  real. 

Having  derived  the  decay  rates  in  parabolic  coordinates, 
we  now  define  a  more  natural  set  of  parameters  to  describe 
the  system  [see  Fig.  3(b)].  Let  us  introduce  a  length  scale  w 
that  characterizes  the  curvature  of  the  nanotip  via  p(z) 
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=  ^fwz  (z  ^  0),  where  p  is  the  radius  of  the  nanotip  at  position 
z  (note  also  the  corresponding  shift  in  the  apex  of  the  tip 
from  z=-Vq/ 2  to  z=0).  Furthermore,  let  z--d <  0  be  the 
position  of  the  emitter  [d=(v,2-v q)/2  is  the  distance  be¬ 
tween  the  emitter  and  end  of  the  nanotip].  In  terms  of  these 
parameters,  the  spontaneous  emission  rates  derived  above  for 
a  dipole  p0  z  positioned  along  the  z  axis  can  be  rewritten  as 


+  [l  +(4  dlw)Y 


(45) 


-Tnonrad  _  5 _ 1  ,  /  ^2  ^1  j 

To  ~^'2(k0dZ  m\e2+eJ’ 


(46) 


Eel  - 

= 


1 


/sTi[cVl  +  (4c//vv)J,  (47) 


To  p'(A:0h’)3[1 +  (4J/w)] 
where  apl  only  depends  on  el2  and  is  given  by 
24tt  Aex- e2)Ix{CMC) 


Ot  «1  = 


63/2 


dx(C)/dx 


(48) 


As  in  the  nanowire  case,  one  can  define  a  Purcell  factor  P 
=  rpl/r'  for  the  nanotip,  which  depends  on  the  curvature 
parameter  w  and  emitter  position  d.  Optimization  of  the  Pur¬ 
cell  factor  over  the  emitter  position  yields  a  large  enhance¬ 
ment  of  P  ~  2.5  X  103  as  w— » 0,  which  is  limited  ultimately 
by  dissipation  in  the  metal.  In  Fig.  4(c),  we  plot  P  as  func¬ 
tions  of  w  and  d.  Once  again,  it  can  be  seen  that  achieving 
large  Purcell  factors  does  not  depend  too  sensitively  on  the 
emitter  position. 

To  compare  the  benefits  of  the  nanotip  over  the  nanowire, 
the  relevant  quantity  is  not  just  the  Purcell  factor  but  how 
efficiently  the  generated  plasmons  can  propagate  from  the 
end  to  some  larger  radius  R  further  down  the  tip  [see  Fig. 
3(b)].  Because  the  plasmon  modes  here  were  obtained 
through  a  quasistatic  approximation,  this  calculation  yields 
no  information  about  dissipative  losses  as  the  plasmon 
propagates  along  the  nanotip.  For  example,  in  this  limit  H 
~  0,  so  one  cannot  obtain  the  Poynting  vector  for  the  system. 
To  estimate  the  effect  of  propagative  losses,  however,  we  can 
make  an  eikonal  approximation,39  assuming  that  the  plas¬ 
mons  are  emitted  completely  into  the  end  of  the  tip  (z=0) 
and  that  the  propagative  losses  thereafter  at  any  position  z 
are  identical  locally  to  those  of  a  nano  wire  of  radius  p(z). 
Specifically,  the  fraction  of  the  total  plasmon  emission  that 
successfully  propagates  from  the  end  to  some  larger  radius  R 
is  estimated  to  be 


F  /  r*(R)  \ 

—  =  exp|  -  2  J  Im  klp(z)]dzj .  (49) 

Whereas  Tpl  is  the  spontaneous  emission  rate  into  the  plas¬ 
mons,  f  vi(R)  is  the  rate  that  these  emitted  plasmons  success¬ 
fully  propagate  to  radius  R.  It  is  also  convenient  to  define  an 
effective  decay  into  other  channels,  f'(/?)  =  r'  +  [rpl 

-Tpl(/?)],  which  includes  the  rate  at  which  emitted  plasmons 
are  dissipated  before  reaching  R.  Similarly,  we  define  an  ef¬ 


fective  Purcell  factor  for  the  nanotip,  given  by  P(R) 

=  fpl(/?)/r '(/?),  which  characterizes  the  efficiency  that  plas¬ 
mons  are  generated  and  propagate  successfully  to  R.  In  gen¬ 
eral,  this  quantity  depends  both  on  the  emitter  position  and 
the  curvature  of  the  nanotip.  In  Fig.  4(a),  we  plot  P(R)  when 
optimized  over  w  and  d.  It  can  be  seen  that  the  effective 
Purcell  factor  for  the  nanotip  compares  favorably  to  the  Pur¬ 
cell  factor  of  a  nano  wire  when  k^R^  0.05.  In  other  words, 
the  nanotip  makes  it  possible  to  use  these  larger  sizes  to 
reduce  dissipation,  without  as  steep  a  trade-off  in  the  effi¬ 
ciency  of  coupling  to  the  plasmons.  We  note  that  the  nanotip 
system  also  has  the  added  benefit  of  generating  guided  plas¬ 
mons  along  a  single  direction  of  propagation  away  from  the 
end  of  the  tip. 

We  now  discuss  the  limits  of  validity  of  the  equations 
derived  above  for  the  nanotip.  First,  as  discussed  in  Sec. 
Ill  A,  the  quasistatic  calculation  of  the  decay  rates  is  valid 
only  when  one  is  considering  length  scales  where  the  elec¬ 
trodynamic  Green’s  function  can  be  approximated  by  its 
electrostatic  equivalent.  For  the  nanotip,  this  restricts  the  re¬ 
gimes  of  validity  to  \k\w,  \ki\d<l.  An  additional  set  of  as¬ 
sumptions  is  made  in  using  the  eikonal  approximation  to 
arrive  at  Eq.  (49).  Specifically,  we  have  assumed  that  the 
change  in  the  radius  p(z)  of  the  nanotip  occurs  slowly 
enough  that  the  wave  vector  k\\(p )  of  the  surface  plasmon 
modes  can  adiabatically  follow  the  nanowire  solution  at  each 
position  z-  The  adiabaticity  condition  can  be  quantified  by  a 
parameter  /3= d[l/Re  k\\(p)]/dz,  which  must  remain  small  at 
all  positions  z.  Assuming,  for  example,  that  we  are  consider¬ 
ing  sufficiently  small  length  scales  that  the  wave  vector 
kn[p(z)]^C/p(z)  follows  quasistatic  behavior,  the  condition 
(3<  1  is  satisfied  only  in  the  region  z  >  zc,  where  zc 
=  w/|C2|  represents  some  crossover  point  between  adiabatic 
and  nonadiabatic  propagations.  Thus,  at  first  glance,  the  ac¬ 
cumulated  losses  predicted  by  Eq.  (49)  in  propagating  from 
the  end  of  the  tip  (z  =  0)  to  z~zc  appear  invalid.  We  now 
argue,  however,  that  as  long  as  the  upper  limit  of  integration 
z(R)  satisfies  z(R)  >  zc,  Eq.  (49)  remains  a  good  approxima¬ 
tion  for  finding  the  total  losses.  In  particular,  we  argue  that 
the  actual  losses  from  z=0  to  z~zc  are  quite  small,  com¬ 
pared  to  the  losses  accumulated  in  going  from  z~zc  to  z 
>  zc.  In  the  latter  region,  the  integrand  of  Eq.  (49)  is  a  good 
approximation,  as  previously  discussed.  Then,  in  the  region 
from  z=0  to  z~zc,  the  specific  form  of  the  integrand  does 
not  matter,  as  long  as  it  provides  a  small  overall  contribution 
to  the  total  losses.  First,  it  is  straightforward  to  show  that  the 
predicted  loss  in  this  region,  while  incorrect,  is  small  for 
realistic  systems,  rpl[7?(zc)]/rpl~exp(-Im  C/|C|).  Next,  we 
argue  that  the  actual  losses  in  this  region  should  be  small. 
Physically,  the  quantity  \qo\~2  =  vl/\C\2~  zc  determines  the 
fundamental  length  scale  for  the  nanotip  plasmons,  since  q0 
is  the  wave  vector  for  the  plasmons  in  parabolic  coordinates. 
Because  dissipation  can  be  treated  as  a  small  perturbation  for 
realistic  systems,  the  losses  should  be  significant  only  on 
length  scales  much  larger  than  this,  just  as  losses  are  only 
relevant  on  scales  much  longer  than  the  inverse  plasmon 
wave  vector,  1  /  |^J,  for  the  nano  wire. 

To  check  the  analytical  results  derived  above  for  the  nan¬ 
otip,  we  have  also  performed  detailed  numerical  simulations 
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FIG.  5.  (Color  online)  Numerically  calculated  fields  due  to  a  dipole  emitter  near  a  conducting  nanotip,  obtained  by  the  boundary  element 
method,  (a)  The  energy  flux  |Re(EXH*)|,  in  arbitrary  units.  The  position  of  the  emitter  is  denoted  by  the  blue  circles,  while  the  boundary 
of  the  nanotip  is  given  by  the  dotted  lines.  The  plots  shown  are  for  a  final  nanotip  radius  of  k0R  =  03,  curvature  parameter  k0w= 0.022,  and 
emitter  positions  k0d= 0.002,  0.2,  and  0.7.  It  can  be  seen  that  both  the  total  spontaneous  emission  rate  Ttotal  and  the  emission  rate  into 
plasmons  increase  as  the  emitter  is  brought  closer  to  the  nanotip.  (b)  The  quantity  |Re(E  X  H*)  |  /rtotal,  for  the  same  parameters.  This  quantity 
is  proportional  to  the  energy  flux  normalized  by  the  total  power  output  of  the  emitter.  The  k$d= 0.002  plot  is  mostly  dark,  indicating  that  most 
of  the  decay  is  into  nonradiative  channels.  The  k0d=0.2  case  is  characterized  by  bright  spots  along  the  entire  edge  of  the  nanotip,  which 
indicates  efficient  plasmon  excitation.  The  k0d= 0.7  case  exhibits  the  typical  lobe  pattern  associated  with  radiative  decay  of  a  dipole. 


using  the  boundary  element  method  (BEM).40  Details  of  our 
implementation  are  given  in  Appendix  D.  BEM  simulations 
are  fully  electrodynamic  solvers  of  Maxwell’s  equations,  and 
they  were  used  to  obtain  the  classical  electromagnetic  field 
solutions  of  an  oscillating  dipole  emitter  p0e~lo)t  near  a  nan¬ 
otip.  The  results  of  a  few  sample  simulations  are  shown  in 
Fig.  5,  for  a  tip  curvature  parameter  k0w= 0.022,  final  radius 
k0R=03,  and  varying  emitter  positions  k0d= 0.002,  0.2,  and 
0.7.  In  Fig.  5(a),  we  plot  the  quantity  |Re(E  X  H*)|,  which  is 
proportional  to  the  Poynting  vector  and  corresponds  to  the 
total  energy  flux  of  the  system.  The  nanotip  is  assumed  to  be 
composed  of  silver  in  a  surrounding  dielectric  e1  =  2,  and  its 
boundary  is  given  by  the  dotted  line,  while  the  emitter  posi¬ 
tions  are  denoted  by  circles.  The  total  spontaneous  emission 
rate  is  given  via  rtotal=(rpl+r')ocIm[p0-E1(r' ,r')]  and  is 
determined  numerically  for  each  configuration  by  finding  the 
total  field  at  the  dipole  location.  On  the  other  hand,  the  ef¬ 
fective  emission  rate  f pl(R)  into  the  plasmons  is  determined 
by  a  best  fit  of  the  fields  in  the  region  of  constant  R  to  the 
known  plasmon  solution  on  a  nano  wire  given  in  Eq.  (3),  and 
then  by  calculating  the  total  power  transport  of  this  best-fit 
mode  through  the  integrated  Poynting  vector.  This  total 

power  is  directly  proportional  to  fpl(/?).  The  figure  confirms 
the  qualitative  behavior  that  we  expect  and  have  described 
previously.  In  particular,  the  generated  plasmon  field  and  to¬ 
tal  spontaneous  emission  rate  are  largest  for  very  small  sepa¬ 
rations  and  decrease  as  the  emitter  is  placed  further  away 


from  the  end  of  the  nanotip.  In  Fig.  5(b),  we  plot  |Re(E 
XH*)|/rtotal,  which  is  proportional  to  the  energy  flux  nor¬ 
malized  by  the  total  power  output  of  the  emitter.  This  quan¬ 
tity  yields  information  about  the  efficiency  of  decay  into  the 
various  channels.  For  small  separations  (k0d= 0.002),  the  plot 
is  mostly  dark,  which  indicates  that  the  decay  of  the  dipole  is 
predominantly  nonradiative.  For  k0d=  0.2,  the  maximum 
(corresponding  to  bright  spots  in  the  plot)  is  located  along 
the  entire  surface  of  the  nanotip,  which  indicates  highly  ef¬ 
ficient  plasmon  excitation.  Here,  although  the  total  emission 
rate  into  plasmons  decreases  from  the  k0d= 0.002  case  [as 
seen  in  Fig.  5(a)],  the  efficiency  increases  dramatically  due 
to  less  competition  from  nonradiative  decay.  Finally,  for 
k0d= 0.7,  the  maximum  appears  as  the  typical  lobe  pattern 
associated  with  radiative  decay. 

In  Fig.  4(a),  we  have  plotted  the  numerically  optimized 

values  of  P(R)  for  a  few  values  of  R.  It  can  be  seen  that  the 
values  obtained  through  analytical  approximations  and  nu¬ 
merical  BEM  closely  agree.  Unlike  the  theoretical  predic¬ 
tions,  however,  the  numerically  calculated  error  probability 
does  not  increase  monotonically  with  R.  We  believe  that  the 
origin  of  this  is  that,  for  the  numerically  optimized  param¬ 
eters,  the  condition  R  >  p(zc)  under  which  the  theoretical  pre¬ 
dictions  hold  is  only  weakly  satisfied,  and  the  excitation  re¬ 
gion  for  the  plasmons  cannot  strictly  be  thought  of  as  a 
single  point  at  the  end  of  the  tip  (z=0). 
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V.  SINGLE-PHOTON  GENERATION  VIA  COUPLING 
TO  DIELECTRIC  WAVEGUIDE 

We  have  shown  in  previous  sections  that  a  single  emitter 
can  spontaneously  emit  into  the  guided  plasmon  modes  of  a 
nearby  nanostructure  with  high  probability.  This  prospect  of 
efficient  conversion  between  an  excitation  of  the  emitter  and 
a  single  photon  has  a  number  of  applications  in  the  fields  of 
quantum  computing  and  quantum  information.  In  this  sec¬ 
tion,  we  consider  one  particular  application  involving  the  use 
of  such  a  system  as  an  efficient  single-photon  source.  The 
concepts  behind  single-photon  generation  on  demand  with 
an  individual  emitter  in  a  cavity  have  been  discussed 
elsewhere41-43  and  will  not  be  presented  in  detail  here.  We 
note  also  that  the  ideas  behind  single-photon  sources  can  be 
extended  to  create  long-distance  entanglement  between  emit¬ 
ters,  as  detailed,  e.g,  in  Ref.  44. 

Because  of  dissipative  losses  in  metals,  the  plasmon 
modes  are  not  directly  suitable  as  carriers  of  information 
over  long  distances.  We  show,  however,  that  plasmonic  de¬ 
vices  can  serve  as  an  effective  intermediate  step,  and  in  par¬ 
ticular,  can  be  efficiently  outcoupled  to  the  modes  of  a  co¬ 
propagating  dielectric  waveguide.  The  single-photon  device 
is  illustrated  schematically  in  Fig.  1.  In  Fig.  1(a),  an  optically 
addressable  emitter  with  multiple  internal  levels  sits  in  the 
vicinity  of  a  conducting  nano  wire.  The  emitter  is  strongly 
coupled  to  the  nanowire,  such  that  single  photons  on  demand 
can  be  generated  with  high  efficiency  in  the  plasmon  modes 
by  external  manipulation  of  the  emitter.  The  addressability  of 
the  emitter  along  with  the  internal  levels  allows  for  shaping 
of  this  single-photon  pulse,45  as  illustrated  in  Fig.  1(b).  Here, 
a  three-level  emitter  is  shown  with  two  ground  or  metastable 
states  Is)  and  | g).  We  assume  that  the  transition  between  | s) 
and  excited  state  \e)  is  decoupled  from  the  plasmon  modes, 
due  to,  e.g,  the  orientation  of  its  dipole  moment,  but  that  the 
states  are  coupled  via  some  classical  control  field  with  Rabi 
frequency  H(r).  We  also  assume  that  the  transition  \e)-\g)  is 
coupled  via  the  plasmon  modes  of  the  nano  wire;  i.e.,  the 
state  \e)  can  decay  at  a  rate  Tpl  into  | g)  by  emitting  a  photon 
into  the  plasmon  modes.  In  addition,  there  is  a  small  rate  T' 
at  which  the  excited  state  can  decay  without  emitting  a  plas¬ 
mon.  A  single  photon  in  the  plasmon  modes  of  the  nanowire 
is  generated  with  high  probability  by  initializing  the  emitter 
in  Is)  and  exciting  the  transition  \s)-\e)  with  the  control  field 
H(^).  The  decay  of  \e)  into  | g)  and  the  generation  of  an 
outgoing  single  plasmon  then  occur  with  high  probability, 
with  the  shape  of  the  single-photon  wave  packet  determined 
by  the  shape  of  H(r).  We  further  assume  that  the  plasmon  is 
then  evanescently  coupled  to  the  nearby  dielectric  waveguide 
[see  Fig.  1(a)],  which  copropagates  with  the  nanowire  over 
some  distance  Lex  during  which  this  coupling  is  non- 
negligible.  The  coupling  is  a  reversible  process,  and  the  dis¬ 
tance  Lex  is  optimized  to  maximize  efficiency  of  ending  up 
with  a  single  photon  in  the  waveguide  (i.e.,  to  prevent  further 
Rabi  oscillations  back  into  the  nano  wire).  A  similar  setup 
with  a  nanotip  is  illustrated  in  Fig.  1(c).  Here,  the  nanotip 
radius  p(z)  expands  to  some  final  radius  R ,  at  which  point 
coupling  with  the  waveguide  starts  to  occur.  Initiating  the 
coupling  once  the  nanotip  has  reached  a  constant  radius  al- 


FIG.  6.  Wave  vector  k\\  of  the  fundamental  guided  modes  of  a 
cylindrical  dielectric  waveguide  with  core  permittivity  ec=  13  and 
surrounding  permittivity  e1=2,  plotted  as  a  function  of  core  radius 


lows  the  two  systems  to  be  easily  mode  matched,  as  dis¬ 
cussed  below.  When  optimized,  we  estimate  that  single¬ 
photon  generation  efficiencies  exceeding  —95%  are  possible 
in  this  tiered  configuration. 

To  treat  the  problem  analytically,  we  consider  the  simple 
situation  of  our  nanostructure  coupled  to  a  cylindrical  dielec¬ 
tric  waveguide  (e.g,  an  optical  fiber)  of  radius  Rg,  such  that 
the  modes  can  be  calculated  analytically  using  the  methods 
described  in  Appendix  A.  It  can  be  shown  that  the  fundamen¬ 
tal  modes  of  the  waveguide  are  degenerate  m-±  1  modes 
that  are  not  cut  off  as  Rg-+0.  The  dependence  of  their  wave 
vector  k\\  on  Rg  is  shown  in  Fig.  6  for  a  core  permittivity 
eg=  13  and  surrounding  permittivity  ex-2.  These  parameters 
correspond  closely  to  that  of  a  Si/Si02  guide  at  X0=l  /zn l. 
To  simplify  the  calculation,  we  also  assume  that  coupling 
between  the  wire  and  higher-order  waveguide  modes  is  neg¬ 
ligible.  This  can  be  achieved,  for  example,  by  operating  be¬ 
low  the  cutoff  radius  of  higher-order  modes  or  by  operating 
with  sufficiently  large  wave  vector  mismatch  between  the 
plasmon  and  higher-order  guide  modes. 

We  make  the  ansatz  that  the  total  field  of  the  system  is 
given  by  a  superposition  of  the  unperturbed  modes  of  the 
nanostructure  and  waveguide.  While  this  cannot  strictly  be 
correct,  as  such  a  solution  violates  boundary  conditions  at 
each  interface,  we  rely  on  such  an  assumption  to  give  us  the 
correct  qualitative  behavior  without  resorting  to  more  com¬ 
plex  numerical  calculations.  Specifically,  we  assume  that  the 
total  electric  field  for  the  system  takes  the  form 

Er(r)=  2  S  (50) 

/ JL=W,g  i=  1 

where  /z  indexes  the  nanostructure  (w)  and  waveguide  (g) 
systems,  and  i=  1 , . . .  runs  over  the  modes  of  system  /z. 
In  the  following,  we  will  explicitly  treat  the  nanotip  case, 
where  the  plasmons  propagate  in  a  single  direction,  although 
this  argument  can  easily  be  extended  to  the  nano  wire.  We 
emphasize  that  we  are  considering  coupling  of  the  plasmon 
mode  to  the  waveguide  once  the  nanotip  has  already  ex- 


035420-12 


STRONG  COUPLING  OF  SINGLE  EMITTERS  TO  SURFACE... 


PHYSICAL  REVIEW  B  76,  035420  (2007) 


panded  to  its  final  radius  R ,  at  which  point  the  plasmon  mode 
solution  becomes  identical  to  that  of  a  nanowire.  In  our  case, 
Nw=  1  as  we  only  consider  the  fundamental  plasmon  mode  of 
the  nanotip,  while  Ng=2  as  we  take  into  account  the  two 
degenerate,  copropagating  fundamental  modes  of  the  wave¬ 
guide.  E^r)  here  represents  the  unperturbed  solution  of 
mode  i  in  system  fi  (without  the  presence  of  the  other  sys¬ 
tem).  A  similar  expression  holds  for  the  total  magnetic  field. 

With  the  ansatz  of  Eq.  (50)  for  the  total  field  of  the  com¬ 
bined  waveguide  and  nanotip  system,  one  can  derive  exact 
equations  of  evolution46  based  on  Lorentz  reciprocity  for  the 
coefficients  C^j.  Explicitly  separating  out  the  plane- wave 
dependence  of  the  unperturbed  fields,  EJbi  i(r)=E/iL  i(p)elk^iZ, 
the  Nw+Ng  coupled-mode  equations  take  the  form 

Nv  dc  _  Nv 

2  ^  P/u,,i-v,M)—M  =  i(oe 2^/w>j(z)C',y(z), 
v=w,g  7=1  dz  v=w,g  7=1 

(51) 

as  derived  in  detail  in  Appendix  E.  The  coefficients  to  the 
system  of  equations  above  are  given  by 

p»,i-,vj(z)  =  ei{k^J~k\w)z  J  dp[EvJ(p)  X  HV(p) 

+  E]jp)  X  HvJ(p)]  ■  l  (52) 

K*i;v/z)  =  J  dpE.Jp)  •  E V(p)[er(p)  -  ev(p) ], 

(53) 

where  eT{p)  is  the  electric  permittivity  of  the  combined  sys¬ 
tem.  Clearly,  the  presence  of  the  phase  factors  el^v’j~k\w^z  in 
the  equations  above  indicates  that,  at  least  under  weak  cou¬ 
pling,  significant  power  transfer  between  the  two  systems 
will  not  take  place  unless  the  two  systems  are  approximately 
mode  matched  with  respect  to  k\\.  In  practice,  this  implies  that 
for  a  final  tip  radius  R ,  there  is  some  ideal  waveguide  size  Rg 
that  allows  for  maximum  transfer  efficiency  between  the  two 
systems.  A  similar  optimization  of  the  waveguide  parameters 
exists  in  the  case  of  arbitrary  coupling  strength  between  the 
two  systems,  although  this  problem  is  more  complex  because 
one  must  account  for  factors  such  as  the  phase  shift  of  one 
system  due  to  the  other.  We  emphasize  that  the  coupled¬ 
mode  equations  above  are  exact  within  the  ansatz  of  Eq. 
(50).  For  example,  for  two  lossless  systems,  these  equations 
conserve  power,  and  for  a  lossy  system  (such  as  a  nanotip), 
the  effects  of  losses  are  treated  exactly.  By  convention,  the 
normalization  of  fields  is  such  that  the  integrals  appearing  in 
the  diagonal  matrix  elements  are  set  to  1. 

For  the  waveguide  and  nanotip  systems  coupled  over  a 
length  Lex ,  the  exact  single-photon  generation  efficiency  will 
depend  on  the  details  of  how  the  two  systems  are  brought 
together  and  separated  apart.  In  practice,  for  example,  the 
two  systems  should  be  brought  together  slowly  enough  that 
the  introduction  of  the  waveguide  does  not  cause  significant 
backscattering  of  the  plasmon,  yet  quickly  enough  that  this 
introduction  length  is  small  compared  to  the  plasmon  decay 
length.  Furthermore,  in  reality,  the  coupling  region  will  not 


be  a  step  of  length  Lex ,  but  will  be  characterized  by  some 
smooth  transition.  To  avoid  the  many  details  associated  with 
this  introduction  and  separation  and  to  approximately  calcu¬ 
late  the  efficiency,  we  will  consider  an  idealized  system  and 
make  three  assumptions  as  follows: 

(i)  The  decay  rates  of  the  emitter  are  not  affected  by  the 
presence  of  the  nearby  dielectric  waveguide.  In  particular, 
the  Purcell  factors  calculated  earlier  for  the  nanotip  are  un¬ 
changed. 

(ii)  The  radius  of  the  nanotip  is  given  by  p(z)  =  Vwz  for 
z<Zo  and  becomes  constant,  R  =  p(z0)  =  \Zwz0,  for  z^Zo-  For 
z^z0,  the  plasmon  mode  solution  becomes  identical  to  the 
nano  wire  solution,  and  in  particular,  has  well-defined  k\\, 
which  allows  it  to  be  easily  mode  matched  with  the  wave¬ 
guide.  It  is  assumed  that  coupling  between  the  nanotip  and 
waveguide  begins  at  z=Zo,  with  the  initial  field  amplitudes  of 
the  coupled  system  given  by 


Cw(zo) 


HR) 

_  l  +  P(R)  _ 


1/2 


(54) 


C*fi(z0)  =  0,  (55) 

where  the  probability  P(7?)/[l  +P(R)~\  of  finding  a  plasmon 
excitation  at  R  is  already  optimized  over  the  nanotip  curva¬ 
ture  and  emitter  position. 

(iii)  Equation  (51)  exactly  describes  the  coupling  between 
the  two  systems  in  the  region  zo^z^Zo+Lex-  To  estimate  the 
probability  of  transfer  from  nanowire  to  waveguide  after  dis¬ 
tance  Lex  when  the  two  systems  are  once  again  separated,  we 
project  the  total  field  of  Eq.  (50)  at  z=Zo+Lex  into  the  wave¬ 
guide  mode.  Specifically,  the  projected  field  amplitude  in  the 
waveguide  in  either  of  the  degenerate  modes  i  is  given  by 

Cproj/z,)  +  Lex)  =  2  J  dp[ Er(r)  X  H*  ,(r)]  •  £,  (56) 

where  the  factor  of  2  arises  due  to  the  normalization  conven¬ 
tion  for  the  unperturbed  modes. 

Because  of  the  symmetry,  the  projected  field  strengths 
|Cproj,;|2  calculated  above  are  equal  for  the  two  degenerate 
waveguide  modes,  and  the  quantity  2|Cproj?/|2  then  corre¬ 
sponds  to  the  efficiency  of  single-photon  generation.  Here, 
the  additional  factor  of  2  accounts  for  the  mode  degeneracy. 
This  quantity  takes  completely  into  account  the  propagative 
losses  of  the  plasmons,  imperfect  coupling  between  the  nan¬ 
otip  and  waveguide,  and  the  Purcell  factor  of  the  nanotip. 

In  Fig.  7(a),  we  plot  the  efficiency  of  single-photon  gen¬ 
eration  as  a  function  of  R  for  both  the  nanowire  and  nanotip 
systems.  For  each  R,  the  plotted  efficiencies  have  been  opti¬ 
mized  over  all  other  possible  parameters  of  the  system.  For 
the  nanowire  configuration,  we  have  assumed  that  the  result¬ 
ing  forward-  and  backward-propagating  waves  in  the  wave¬ 
guide  can  be  perfectly  combined.  In  the  figure,  we  have  also 
included  points  obtained  by  our  BEM  simulations  of  a  nan¬ 
otip.  Here,  we  have  taken  the  numerically  optimized  values 

of  P(R)  and  plugged  them  in  as  initial  values  for  the 
coupled-mode  theory  above.  It  can  be  seen  that  the  numeri¬ 
cal  simulations  agree  well  with  our  theoretical  predictions. 
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FIG.  7.  (a)  Optimized  efficiencies  of  single-photon  generation 
vs  R.  We  have  assumed  that  coupling  to  waveguide  modes  other 
than  the  fundamental  mode  is  negligible,  i.e.,  the  waveguide  is  ef¬ 
fectively  in  the  single-mode  regime.  Solid  (dashed)  line:  theoretical 
efficiency  using  a  nanowire,  using  spontaneous  emission  rates  ob¬ 
tained  by  quasistatic  (fully  electrodynamic)  calculations.  Dotted 
line:  theoretical  efficiency  using  a  nanotip.  Solid  points:  nanotip 
efficiency  based  on  boundary  element  method  simulations,  com¬ 
bined  with  coupled-mode  equations,  (b)  Optimal  coupling  length 
Lex  for  a  nanotip  as  a  function  of  R.  Here,  Lex  is  given  in  units  of  the 
plasmon  wavelength  Xpl  at  that  particular  R. 

We  find  that  photon  efficiencies  of  nearly  80%  are  possible 
for  the  nano  wire, 37  while  efficiencies  exceeding  95%  are 
possible  for  the  nanotip.  In  Fig.  7(b),  we  plot  the  optimal 
coupling  length  Lex ,  in  units  of  Xpl,  as  a  function  of  R  for  the 
nanotip  (. Lex  for  the  nanowire  should  be  twice  that  of  the 
nanotip  to  account  for  the  transfer  of  the  forward-  and 
backward-propagating  components  of  the  emitted  plasmon). 
It  can  be  seen  that  the  outcoupling  to  the  waveguide  can,  in 
principle,  occur  quite  rapidly,  over  length  scales  of  a  few  Xpl. 

The  existence  of  an  optimum  R  for  photon  generation  can 
be  intuitively  understood.  For  smaller  R ,  the  coupling  be¬ 
tween  the  emitter  and  plasmon  modes  can  be  quite  large. 
However,  these  tightly  confined  plasmon  modes  are  accom¬ 
panied  by  higher  propagative  losses  and  cannot  be  as  effi¬ 
ciently  coupled  to  the  waveguide  system.  The  coupling  effi¬ 
ciency  between  plasmons  and  waveguide  modes  improves 
for  larger  R.  For  the  nano  wire,  however,  the  larger  radius 
results  in  weaker  coupling  between  the  plasmon  and  emitter, 
while  for  the  nanotip,  the  accumulated  propagative  loss  in¬ 
creases  as  the  final  radius  grows. 


VI.  EFFECTS  OF  SURFACE  ROUGHNESS 

In  previous  sections,  we  have  treated  the  problem  of  plas¬ 
mon  propagation  on  smooth  nanowires  and  nanotips,  taking 
into  account  inherent  dissipative  losses  characterized  by 
Im  62.  In  practice,  however,  these  structures  are  not  perfectly 
smooth  and  the  surface  roughness  can  give  rise  to  new  scat¬ 
tering  mechanisms  for  the  plasmons.  While  the  general  solu¬ 
tion  for  the  fields  in  the  presence  of  arbitrary  roughness  is  a 
complicated  problem,  we  calculate  the  effects  in  two  limits. 
In  Sec.  VI  A,  we  calculate  the  losses  on  a  nanowire  due  to 
radiative  scattering  in  the  limit  of  small  roughness  and  zero 
heating  (Im  62  =  0).  Here,  the  plasmons  experience  no  inher¬ 
ent  loss  due  to  the  metal  but  can  receive  momentum  kicks 
from  the  roughness  that  cause  them  to  scatter  radiatively.  In 


Sec.  VI B,  we  calculate  the  effects  of  small  roughness  for  a 
nanowire  in  the  nonretarded  limit,  where  radiative  effects  are 
ignored  but  the  effects  of  increased  dissipative  losses  are 
treated.  While  we  explicitly  treat  only  the  nanowire  case 
here,  we  note  that  the  results  obtained  can  also  be  incorpo¬ 
rated  into  our  model  for  nanotip  losses  via  Eq.  (49). 


A.  Radiative  losses 


For  simplicity,  we  consider  a  wire  with  axial  symmetry, 
but  with  a  surface  profile  given  by  p0(z)=R+p£(z),  where  R 
is  the  average  radius  of  the  wire,  £(z)  is  some  random  func¬ 
tion  describing  the  roughness,  and  p  is  an  expansion  param¬ 
eter  that  will  be  taken  to  equal  1  at  the  end.  We  will  calculate 
in  perturbation  theory  the  radiated  field  scattered  from  the 
roughness,  from  an  initial  field  corresponding  to  the  funda¬ 
mental  plasmon  mode  for  a  perfectly  smooth  wire.  Because 
of  the  symmetry,  the  only  nonzero  components  of  the  fields 
remain,  Ep,  Ev  and  H^,  which  will  also  have  axial  symmetry. 
As  will  be  seen  later,  it  suffices  for  now  to  consider  only  Ez, 
as  the  other  components  depend  on  Ez  in  a  simple  way 
through  Maxwell’s  equations.  We  proceed  by  breaking  up  the 
total  field  along  z  in  region  i  into  incident  and  scattered  fields 


rr total  _  /7°  _i_  fs 

ni,z  ~  ^ i,z  +  ^i,z’ 


(57) 


where  E®z  is  the  z  component  of  the  fundamental  plasmon 
mode  given  by  Eqs.  (3)  and  (4),  and  further  assume  that  the 
scattered  field  can  be  expanded  in  a  power  series 


^=2^.  (58) 

n=  1 

In  the  following,  we  will  calculate  the  first-order  scattered 
field  E^.  We  make  the  ansatz  that  E^  can  be  expanded  in 
the  form47 

f00  h 2 

£(i!L  dh^ih^p^Aih^, 

J  —oo  k] 

r  h2 

4!z=  ^Vo^2±p)-^(6||VY  (59) 

J  —00  1^2 

where  each  Fourier  component  is  an  outgoing  solution  of  the 
wave  equation  with  appropriate  boundary  conditions  at  p 
=  0  and  p=°°,  as  derived  in  Eq.  (A4).  From  Eq.  (A4),  one 
also  sees  that  Ep  and  H ^  are  determined  completely  once  Ez 
is  known.  Using  these  relations,  the  total  (incident  plus  scat¬ 
tered)  fields  Etotal  and  Htotal  are  straightforward  but  lengthy 
to  write  down,  and  are  given  to  order  p  in  Eq.  (FI)  in  Ap¬ 
pendix  F. 

The  coefficients  A(/ty)  and  B(h\ |)  are  determined  by  en¬ 
forcing  continuity  of  the  tangential  fields  at  the  boundary 
p0(z).  Specifically,  we  require  that 

(?-Er)u+^=  (t-*rx=R+H(z)’ 


t  = 


z 


+  P~TP 
dz 


1  +p2 


d£ 

dz 


2’ 
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rrtOtal  I  _  rrtOtal  I  ( 

I p=R+p£(z)  -  M<f>,2  I P=R+P&z)> 

where  t(z )  is  the  unit  tangent  vector  to  the  interface.  These 
equations  can  be  solved  perturbatively  by  expanding  them  in 
p  and  solving  each  order.  The  O(p0 )  equation  is  trivially 
satisfied  by  the  fundamental  plasmon  mode  for  a  smooth 
nano  wire.  To  solve  to  0(p ),  it  is  useful  to  first  introduce  the 
Fourier  transform  of  the  surface  roughness, 

£(z)=  [  ? eihH{h{ |).  (61) 

J-  OO  27 7 

Using  the  Fourier  transform  £(Ay),  the  0(p)  equations  be¬ 
come  algebraic  in  Fourier  space  and  have  solutions  (see  Ap¬ 
pendix  F) 


A(A„)  = 


StMiL  ) 

2tt  hi± 


physically  corresponds  to  the  power  radiated  away  to  infin¬ 
ity.  Specifically,  the  radiated  power  per  unit  area  is  given  by 


sp=-1-{e^;)  oo) 

1  r  J(.  „,ihlL2h[\ 
2  co/z0J_ 


(66a) 


dhudh. 


k] 


X{A{h^A\ti«))ei(hrh>. 


H0(hl±p)H^(h[±p) 


(66b) 


Substituting  the  solution  for  A  (Ay)  derived  in  Eq.  (62)  and 
using  the  correlations  in  Eq.  (64),  it  is  straightforward  to 
evaluate  the  integral  over  A, \  and  arrive  at 


^Va3  f‘  dh{le-^a^-k^hn 
4  Vtt-  J-k. 


XH0(hllP)H '  (hllP)\f(h^\ 


(67) 


(62) 

277  ^2_L 

where  k\\  denotes  the  unperturbed  plasmon  wave  vector  (in 
this  section,  we  take  Im  e2  =  0  so  that  k\\  is  purely  real).  The 
scattering  coefficients  /(Ay)  and  g(h^  are  complicated  func¬ 
tions  of  A,,  and  R ,  and  are  given  in  Appendix  F.  Physically, 
the  equations  above  state  that,  to  first  order,  the  surface 
roughness  contributes  single  momentum  kicks  to  the  unper¬ 
turbed  plasmon  fields  with  a  strength  determined  by  the  Fou¬ 
rier  components  of  the  roughness.  From  this  point  forward, 
we  set  p  =  1 . 

We  now  consider  some  random  surface  profile  such  that 

<f(z)>  =  0, 

<£(£)£&))  (63) 

with  corresponding  correlations 

<£(*)>  =  0, 

<; l{k)T(k '))  =  2tj 3l2S2ae-{1,4)a2k2S(k-k'),  (64) 

for  the  Fourier  components.  Physically,  8  and  a  correspond 
respectively  to  the  typical  amplitude  and  length  of  a  rough 
patch  on  the  surface  of  the  wire.  It  is  also  useful  to  define 
5=  8/ a  as  a  typical  “slope”  to  the  roughness.  To  calculate  the 
power  radiated  due  to  the  surface  roughness,  we  find  the 
ensemble-averaged  Poynting  vector  far  from  the  wire.  It  is 
sufficient  to  consider  just  the  component  of  (S)  oriented 
along  p,  given  outside  the  wire  by 

Sp  =  -X-{E«fH*« ?),  (65) 

where  the  fields  and  are  given  to  first  order  by  Eq. 
(FI).  The  calculation  of  Sp  simplifies  further  because  the 
incident  plasmon  field  decays  exponentially  away  from  the 
wire,  and  thus,  to  lowest  order  only  the  first-order  scattered 
fields  will  contribute  to  the  Poynting  vector  at  large  p,  which 


In  the  expression  above,  we  have  truncated  the  bounds  of  the 
integral  to  ±kx  because  we  are  interested  in  the  Poynting 
vector  far  away  from  the  wire,  where  only  radiative  fields 
| A|||  k\  contribute.  With  knowledge  of  the  Poynting  vector, 
it  is  then  possible  to  find  the  dissipation  rate  of  the  plasmons 
due  to  radiative  scattering,  given  by 

r  _  r  _ pjySp _ 

1  rad, rough  \  C  ^ 

~  J  dpe0—[e(p,(o)(o]\E(p)\2  +  p,0\H(p)\2 

(68) 

The  denominator  on  the  right-hand  side  of  the  equation 
above  can  be  identified  with  the  plasmon  energy  per  unit 
length. 

We  first  qualitatively  discuss  the  behavior  of  Trad  rough  be¬ 
fore  deriving  various  limits  more  quantitatively.  From  Eq. 
(67),  it  is  clear  that  rrad  rough  scales  explicitly  like  8 2  or  s2. 
Physically,  this  occurs  because  the  lowest-order  contribution 
to  the  Poynting  vector  far  away  from  the  wire  is  due  to  the 
combination  of  a  first-order  scattered  electric  field  and  first- 
order  scattered  magnetic  field.  In  Fig.  8,  the  quantity 
rrad, rough/ “S2 w  is  evaluated  numerically  as  a  function  of  wire 
radius  R  and  correlation  length  alR  for  a  silver  nano  wire  at 
\0=1  pm  and  e1  =  2.  We  are  particularly  interested  in  the 
nano  wire  limit,  when  the  plasmon  wave  vector  &u~  C/R.  We 
see  that  for  fixed  R ,  the  scattering  reaches  a  peak  for  some 
particular  value  of  a/ R.  More  careful  inspection  reveals  that 
the  maximum  occurs  when  a^R/ C^X pl.  This  result  makes 
intuitive  sense,  since  the  characteristic  momentum  kick 
~  1  /a  that  the  plasmon  wave  vector  k\\  receives  due  to  rough¬ 
ness  must  be  on  the  order  of  Cl R  in  order  for  the  resulting 
wave  vector  to  lie  in  the  radiative  range  between  -kx  and  kx. 
In  the  limit  a/R>C ,  one  observes  an  exponential  suppres¬ 
sion  of  scattering,  due  to  the  fact  that  the  roughness  has  a 
very  narrow  momentum  distribution  and  cannot  possibly 
contribute  a  large  kick  to  k\\.  In  fact,  in  this  regime,  one 
physically  expects  the  plasmon  wave  vector  to  adiabatically 
vary  with  the  changing  wire  radius.  In  the  other  limit  a/R 
<C,  the  scattering  also  decreases,  but  with  a  polynomial 
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FIG.  8.  (Color  online)  The  plasmon  dissipation  rate  due  to  ra¬ 
diative  scattering  off  of  surface  roughness,  rrad  r0Ugh/s2u>,  as  func¬ 
tions  of  wire  radius  R  and  correlation  length  a/ R.  The  numbers  are 
calculated  for  a  silver  nanowire  at  X0=l  jam  and  e1  =  2. 


dependence  on  R,  as  will  be  proven  below.  Here,  the  momen¬ 
tum  distribution  of  the  roughness  becomes  very  wide,  and 
thus,  the  probability  of  receiving  a  kick  that  results  in  a  final 
momentum  between  ±kl  becomes  quite  small.  Finally,  for 
fixed  slope  v,  it  can  be  seen  that  the  scattering  decreases  as 
R^O  at  any  correlation  length  a.  This  result  is  also  easily 
understood,  as  the  plasmon  wave  vector  k\\  becomes  increas¬ 
ingly  far  removed  from  the  range  of  radiative  wave  vectors. 
In  Table  I,  we  calculate  the  scattering  rates  for  wire  sizes 
k0R=OA,  0.2,  and  0.3  (or  16,  32,  and  48  nm),  for  a  few 
chosen  roughness  parameters.  The  scattering  rates  are  given 
as  a  percentage  increase  in  Im  k\\  over  the  values  for  a  smooth 
nanowire.  It  can  be  seen  that  strong  suppression  of  radiative 
scattering  occurs  both  for  smaller  R  and  when  a  is  either 
much  larger  or  much  smaller  than  R ,  which  confirms  our 
earlier  observations.  Furthermore,  it  is  evident  that  under 
reasonable  parameters,  the  losses  in  the  system  are  increased 
only  slightly  due  to  radiative  scattering. 

We  now  analyze  more  carefully  the  behavior  of  the  radia¬ 
tive  scattering  in  the  nano  wire  regime.  To  simplify  the  ex¬ 
pression  further,  we  first  note  that  since  we  are  interested 
in  the  far  field  (p— >o°),  we  can  take  the  asymptotic  limits 
of  the  Hankel  functions  in  Eq.  (67),  H0(h11p)Hl)*(hl_Lp) 


~-2i/(i Thi±p).  One  can  also  derive  an  asymptotic  relation¬ 
ship  of  /(/t||)  as  /?— >0  (see  Appendix  F),  which  upon  substi¬ 
tution  yields 

r rad, rough  =  oe^a3  f  dh„h2xe~(m)a2(kii  -  ft»)2 

V  J-ki 

(*-  0),  (69) 


^^^y^-JodO-Kac).  (70) 

Jo(lC) 

Here,  V  is  a  complicated  dimensionless  parameter,  but  most 
importantly,  only  depends  on  e12- 

From  the  equation  above,  it  is  clear  that  there  are  three 
distinct  regimes  of  interest  defined  by  the  quantity  a  =  k\\a 
=  2i Ta/\pi~  Ca/R,  which  characterizes  the  typical  extent  of 
a  rough  patch  compared  to  the  plasmon  wavelength.  In  the 
limit  a<  1,  one  can  approximate  the  exponential  in  the  inte¬ 
grand  of  Eq.  (69)  as  a  constant,  which  leads  to  straightfor¬ 
ward  evaluation  of  the  integral, 

rrad,rough=  ^ j  O?  (<*«§!).  (71) 


Here,  the  noise  spectrum  of  Eq.  (64)  becomes  very  wide  and 
leads  to  an  a3  scaling  of  the  dissipation  rate.  In  the  opposite 
limit  a>  1,  the  value  of  the  exponential  term  becomes  expo¬ 
nentially  small,  with  a  corresponding  exponential  suppres¬ 
sion  of  the  scattering  rate.  A  more  careful  evaluation  of  the 
integrand  yields 


-  rad, rough  ' 


v 


8tt 

Ca 


■(l/4)a2(£||  -  kx)2 


1) . 


(72) 


Finally,  one  can  show  that  for  fixed,  sub  wavelength  R ,  the 
radiative  scattering  is  most  significant  when  a~0(  1).  In  this 
case,  the  exponential  appearing  in  Eq.  (69)  is  neither  expo¬ 
nentially  small  nor  constant.  However,  one  can  make  the 
rough  approximation  e-(l/4)a2(k\rh\\)2 ^  \ -(l/4)a2(^||-/Z||)2  to 
get  an  idea  of  the  scaling  in  this  regime.  It  is  straightforward 
to  show  that  the  scattering  rate  has  a  maximum  with  respect 


TABLE  I.  Losses  due  to  radiative  scattering  off  of  surface  roughness  for  nanowires  of  varying  sizes  and 
roughness  parameters.  The  scattering  rates  are  given  in  terms  of  the  percentage  increase  in  lmk\\  that  one 
would  expect  over  the  values  for  a  smooth  nano  wire. 


Roughness  parameters 

k0R= 0.1 
{R  ~  16  nm) 

k0R= 0.2 
(R  ~  32  nm) 

koR= 0.3 
(1?  =  48  nm) 

a=0.lR,  <5=0.051?  (5=0.5) 

0.09% 

0.5% 

1.4% 

<7=0.11?,  <5=0.11?  (s=l) 

0.4% 

1.9% 

5.6% 

a=R,  <5=0.051?  (5=0.05) 

0.9% 

4.5% 

12% 

a=5R,  <5=0.051?  (5=0.01) 

2.8% 

8.0% 

10% 

a=l0R,  <5=0.11?  (5=0.01) 

7.0% 

14% 

16% 

a =201?,  <5=0.11?  (5=0.005) 

0.9% 

2.9% 

3.8% 

a =251?,  <5=0.11?  (5=0.004) 

0.3% 

1.3% 

1.8% 
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to  a  at  (12/5)172,  with  a  corresponding  maximum  decay 
rate 


max{rrad ,rough}«  ^Wf]/2J2(^|  .  (73) 

a  y  \  C  / 

Again,  the  radiative  scattering  is  most  significant  when  the 
length  scale  a  of  the  roughness  is  on  the  order  of  the  plasmon 
wavelength,  and  the  maximum  scattering  (for  fixed  s)  de¬ 
creases  as  R  —>  0  due  to  the  increasing  mismatch  between  k\\ 
and  radiative  wave  vectors. 

We  now  consider  the  limits  of  validity  of  the  derivations 
above,  specifically  considering  the  expansions  made  in  Eq. 
(F2)  that  are  necessary  for  the  perturbative  method  used  here. 
The  first  of  these  expansions  requires  that  \d£/dz\  ^  1,  which 
can  be  rewritten  as  a  condition  on  the  slope,  s<l.  Physi¬ 
cally,  this  requirement  states  that  the  typical  length  of  a 
rough  patch  be  much  larger  than  its  typical  height.  The  sec¬ 
ond  line  of  Eq.  (F2)  requires  that  \ki±£\  <  1.  In  the  nano  wire 
regime,  this  requirement  is  equivalent  to  S<\ph  which  states 
that  the  height  of  a  rough  patch  must  be  much  smaller  than 
the  plasmon  wavelength.  Finally,  the  third  line  of  Eq.  (F2) 
requires  \hi±£\<l9  within  the  range  of  hi±  that  are  appre¬ 
ciably  scattered  into.  From  Eqs.  (62)  and  (64),  we  see  that 
the  relevant  range  for  the  parallel  component  of  the  wave 
vector  is  given  by  k\\-  l/a^h\\^k\\  + 1 /a,  and  thus,  the  larg¬ 
est  relevant  transverse  wave  vector  is  K±,maxl 
~max_1^6,^1|^6;(<u/c)2-(&||+  0/a)2 1.  In  the  nano  wire  re¬ 
gime,  k\\~C/R,  there  are  two  limiting  cases.  The  first  is 
when  the  correlation  length  a  is  much  larger  than  R,  a  in 
which  case  \hj  ,  ;Tnax|  —  C/R  and  reduces  to  8<\ph 

In  the  other  limiting  case,  R>a ,  one  finds  that  \hi±^mSLX\ 
~  1  /a  and  the  corresponding  requirement  is  given  by  1. 

We  finally  note  that  while  the  radiative  scattering  goes 
like  S2  or  s2,  the  relevant  quantity  for  dissipative  (heating) 
losses  due  to  roughness  becomes  S  inside  the  wire.  For  this 
quantity,  the  lowest-order  correction  to  the  smooth  wire  so¬ 
lution  will  come  from  a  combination  of  a  first-order  and 
zeroth-order  field.  Thus  one  expects  roughness-induced  dis¬ 
sipative  losses  to  contribute  a  decay  term  proportional  to  8  or 
s,  which  for  small  roughness  will  dominate  over  radiative 
scattering.  This  correction  will  be  treated  in  the  next  section. 


B.  Nonradiative  losses 


satisfy  the  necessary  boundary  conditions  in  the  presence  of 
surface  roughness.  We  can  once  again  find  the  positions  and 
widths  of  the  poles  associated  with  the  system,  which  are 
now  altered  by  the  roughness. 

We  first  write  down  appropriate  expansions  for 
O0r2(r,r').  The  expansion  for  the  incident  component  O0, 
given  originally  in  Eq.  (11),  is  slightly  rewritten  here, 


O0(r,r')  = 


__L _ l_ 

477606!  |r-r' 

00 

1 


:  t~2 - 2  (2-  <5m,o)cosM0-  </>')] 

Zir  606]  m=0 

rco 

x|  dhcos[h(z-z')]Km(hp')Im(hp)  (p<p') 

7o 

:  —2 -  X  e''m(0-0,)  f  dheih{z~z,) 

477  Z6q6j  m—_oo  J—  oo 


XKm(hp')Im(hp)  (p  <  p') . 

(74) 

The  functions  Km(x),Im{x)  are  defined  by 

(75) 

We  also  break  up  Or  2  into  Fourier  components  that  satisfy 
Laplace’s  equation  and  assume  that  these  expressions  hold 
up  to  the  interface: 


<Dr(r,r')  =  -2—  X  eim^  [  dheih^Km(hp)am(h), 
4  7i  6q  6 1  m-_oo  J  —00 

(76) 

1  r°° 

$2(r,r')  =  —j -  X  eim <*“*'>  dheih^Im(hp)/3m(h). 

4  7 T  6q  6 1  m— _oo  J  —00 

(77) 

To  describe  the  surface  roughness,  we  assume  an  interface 
with  axial  symmetry  as  before,  p0(z)=R+p£(z),  where  the 
roughness  profile  £  satisfies  the  correlations  given  in  Eqs. 
(63)  and  (64).  The  coefficients  am  and  (3m  are  determined  by 
the  boundary  conditions,  namely,  that  and  D±  must  be 
continuous  at  the  interface: 


To  study  the  effects  of  surface  roughness  on  nonradiative 
losses,  we  will  make  one  simplifying  assumption  and  calcu¬ 
late  these  losses  in  the  quasistatic  limit.  To  do  this,  we  will 
proceed  in  a  manner  similar  to  that  in  Sec.  Ill  A,  where  we 
found  the  quasistatic  fields  associated  with  a  smooth  nano¬ 
wire.  Here,  the  calculations  for  the  fields  yielded  the  pres¬ 
ence  of  poles  whose  positions  and  widths  give  the  real  and 
imaginary  parts  of  the  wave  vector  k\\.  The  case  of  a  smooth 
nanowire  was  particularly  easy  to  treat  because  of  the  trans¬ 
lational  symmetry  of  the  system.  A  system  containing  surface 
roughness  lacks  such  translational  symmetry,  and  therefore, 
must  be  considered  more  carefully,  but  the  calculation  pro¬ 
ceeds  in  much  the  same  way.  In  particular,  we  will  find  ex¬ 
pressions  for  the  pseudopotentials  O1  =  O0+Or  and  d>2  that 


(**>**  )l p=R+p£(z)  )| p=R+pg(z)’ 

e\n  ■  V$|(r,r')|p=fi+K(s)  =  e2n  ■  V<J>2(r,r')|p=R+pffe), 


n  = 


p-p—z 

dz 


(78) 


Plugging  in  the  expressions  for  O*  given  by  Eqs.  (74),  (76), 
and  (77),  one  can  expand  both  boundary  condition  equations 
to  0(p2).  Then,  replacing  £(z)  in  these  equations  with  its 
Fourier  transform  given  in  Eq.  (61),  one  obtains  a  set  of 
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coupled  equations  for  am  and  /3m(h)  completely  in  Fourier 
space,  which,  unlike  the  case  of  a  smooth  nanowire,  is  not 
decoupled  in  h.  It  is  tedious  but  straightforward  to  show  that, 
to  0(p 2),  this  set  of  equations  is  given  by  the  matrix  integral 
equation 


M0(h) 


+  P 


+  P 


PdJi) 

dqdq' 


I 


da  ~ 

^-ah~q)Ml(h,q) 

2tt 
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(2tt f 
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dqdq' 


£(h-q-  q')C(q')\2(h,q,q')  +  0(p3). 


(79) 


The  matrices  Mt  and  vectors  v*  are  complicated  expressions 
and  are  given  in  Appendix  G.  We  note,  however,  that  in  the 

case  of  no  surface  roughness  (£=0),  the  solution  to  the  re¬ 
sulting  equation  M0(h)  -[am(h)  f3m(h)]T=Vo(h)  reduces  to  that 
of  a  smooth  nano  wire. 

We  now  discuss  how  to  solve  Eq.  (79)  in  the  presence  of 
surface  roughness,  using  the  methods  detailed  in  Ref.  48. 
One  might  first  consider  expanding  am  and  (3m  in  a  power 
series  of  p,  in  a  manner  similar  to  the  field  expansion  in  Eq. 
(58)  for  the  case  of  radiative  scattering,  and  then  solving  the 
0(pn+l)  equations  based  on  the  0{pn)  solutions.  However, 
such  a  perturbative  solution  simply  yields  poles  for  each 
higher-order  correction  with  the  same  location  as  that  of  the 
unperturbed  solutions  and  Mathematically,  this  oc¬ 
curs  because  each  calculation  of  the  next  correction  involves 
an  inversion  On  the  other  hand,  physically,  we  ex¬ 

pect  for  the  surface  roughness  to  result  in  some  shift  of  the 
pole  that  is  not  predicted  by  such  a  perturbative  method.48 
We  thus  consider  an  alternate  approach,  in  which  we  sym¬ 
bolically  sum  the  perturbation  series  in  Eq.  (79)  to  all  orders 
and  then  only  keep  the  lo west-order  result  in  p.  Let  us  sym¬ 
bolically  write  Eq.  (79)  in  the  form 

(Ad0  +  SM)x  =  t0+  St,  (80) 


where  Ad0  and  f0  are  nonrandom  matrices  and  vectors,  re¬ 
spectively,  SM  is  a  random  2X2  matrix  integral  operator,  ST 
is  a  random  vector,  and  x  is  a  column  vector  with  compo¬ 
nents  am  and  f3m.  We  now  define  the  averaging  operator 

Px  =  (x),  (81) 

and  the  operator  Q=  1  -P.  We  can  apply  P  and  Q  to  Eq.  (80) 
to  get 

PM0x  +  PSMx  =  P(f0  +  St),  (82) 


QM0x  +  QSMx  =  Q(  f0  +  St) ,  (83) 

which  after  some  manipulation  results  in  the  set  of  equations 


(A40  +  PSM)(x)  +  PSMQx  =  P(  f0  +  St),  (84) 

Qx  =  (1  +  MQ1QSM)~1Mo1Q(f0  +  St) 

-(1  +MolQSM)-1MolQSM(x).  (85) 

One  can  then  substitute  Eq.  (85)  into  Eq.  (84)  and  solve  for 
(x),  in  which  case  one  obtains 

(A^0  +  ((1  +  SAAMq1  Q)~l  8M)){x) 

=  ((1  +  SMMQ1Q)~l(f0  +  Si)).  (86) 

We  note  that,  unlike  a  perturbative  expansion  and  solution 
for  am  and  (3m,  the  equation  above  is,  thus  far,  exact.  Now, 
we  assume  that  SM  and  St  can  be  expanded  in  powers  of  p 
in  the  form 

SM  —  p  SMI  j  +  p^  SMi  2  +  ...  , 

St  =  pStl+p1St2+  ....  (87) 

Comparing  Eqs.  (79)  and  (87),  we  see  that  (SM\)={Stx) 
=  0  since  (£)  =  0.  With  this  result,  and  utilizing  the  definition 
Q=l-P,  one  can  proceed  to  expand  Eq.  (86)  up  to  0(p2), 
which  yields  (after  setting  p=  1) 

(Ad0  +  (SM2)  -  (SMxMq1  SM\))(x) 

=  (f0)  +  (^2)-(^Ad1Ado1^f1).  (88) 

Substituting  the  corresponding  terms  of  Eq.  (79)  into  the 
equation  above  and  using  the  second-order  correlations  given 
by  Eq.  (64),  we  find  after  simplifying  that 

|m0(/i)  +  J  dq[e~a  q  l4M2(h,h,q) 

_  e-a\h  -  J^“m^  ^ 

=  jv0(/?)  +  ;^=  J  dq[e~a2q2,4\2(h,h,q) 

-  e-a2{h  ~  J ,  (89) 

where  s=S/a. 

We  now  discuss  the  solution  to  a0(h),  which  contains  a 
pole  corresponding  to  the  fundamental  plasmon  mode  m  =  0. 
When  e2  is  a  negative  real  number,  a0  has  a  pole  on  the  real 
h  axis,  whose  position  gives  the  new,  shifted  plasmon  wave 

vector  k\\.  When  e2  has  a  nonzero  imaginary  component,  cr0 
will  have  a  resonance  feature  along  this  axis,  whose  peak 

corresponds  to  Re  k\\  and  whose  width  corresponds  to  Im  k\\. 
A  quick  inspection  of  the  equation  above  reveals  that  k^R 
=  C(€i,s,a/R )  is  a  constant  that  depends  only  on  the  quanti¬ 
ties  en  s ,  and  alR.  Unfortunately,  because  of  the  complexity 
of  Eq.  (89),  it  is  difficult  to  derive  other  scaling  results  for  k\\ 
even  in  limiting  cases.  However,  Eq.  (89)  can  be  solved  nu¬ 
merically.  In  practice,  for  known  parameters,  the  matrices  M 
and  vectors  v  can  be  readily  evaluated  over  some  range  of  h, 
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TABLE  II.  Losses  and  wave  vector  shifts  due  to  nonradiative 
scattering  off  of  surface  roughness  for  nano  wires  with  varying 
roughness  parameters.  The  shifts  in  Re  C  and  changes  in  loss  pa¬ 
rameters  Im  C/Re  C  are  given  in  terms  of  percentage  increase  over 
the  corresponding  values  for  a  smooth  nanowire. 


Roughness  parameters 

A(Re  C) 

A(Im  C/Re  C) 

a=0.1R,  <5=0.011?  (5=0.1) 

0.2% 

0.2% 

a=0.lR,  S=0.05R  (s= 0.5) 

7.5% 

6.8% 

a=R,  ,5=0.011?  (5=0.01) 

0.03% 

1.0% 

a=R,  ,5=0.051?  (5=0.05) 

0.9% 

26% 

a=R,  ,5=0.11?  (5=0.1) 

3.5% 

110% 

a=10R,  <5=0.011?  (5=0.001) 

>0.01% 

2.7% 

a=10R,  <5=0.051?  (5=0.005) 

0.2% 

67% 

a=10R,  <5=0.11?  (5=0.01) 

0.8% 

270% 

from  which  the  solutions  to  the  systems  ct0  and  (30(h)  in  that 
range  immediately  follow.  The  resulting  resonance  in 
Im  a0(h)  as  a  function  of  h  is  then  fitted  to  a  Lorentzian,  with 

its  peak  giving  the  shifted  wave  vector  Re^|  and  its  half¬ 
width  giving  Im  k\\.  In  Table  II,  we  give  the  resulting  losses 
and  wave  vector  shifts  for  a  few  roughness  parameters,  as 
calculated  through  Eq.  (89).  Again  the  numbers  that  we  have 
used  are  for  a  silver  nano  wire  at  \0=1  jmn l  and  ex  =  2.  The 
shifts  in  Re  k\\  (or  equivalently,  Re  C)  and  increases  in  the 
loss  parameter  Im  fcj|/Re  k\\  (or  Im  C/Re  C)  are  given  in  terms 
of  the  percentage  increase  over  their  values  for  a  smooth 
nanowire.  Again  it  can  be  seen  that  for  reasonable  param¬ 
eters,  surface  roughness  adds  only  a  moderate  amount  of  loss 
to  the  system. 


VII.  CONCLUSIONS  AND  OUTLOOK 


nanoscale  traps.  These  traps  might  in  part  be  formed  by  the 
plasmon  fields  themselves. 

We  emphasize  that  the  physical  mechanisms  that  lead  to 
strong  coupling  are  not  restricted  to  the  nanostructures  con¬ 
sidered  here  but  can  be  quite  a  general  feature  of  the  plasmon 
modes  associated  with  subwavelength  conducting  devices.  It 
is  thus  likely  that  the  efficiencies  calculated  here  are  not 
fundamentally  limited  but  can  be  further  improved  by  proper 
design.  Photonic  crystal-like  structures  for  plasmons,54  for 
example,  may  be  a  promising  approach  to  achieve  tight  con¬ 
finement  while  simultaneously  reducing  losses.  Similar 
schemes  may  also  help  to  improve  coupling  between  the 
plasmons  and  dielectric  waveguide  modes.  Such  approaches 
are  likely  to  improve  the  performance  of  plasmon  cavity 
QED  even  further. 
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APPENDIX  A:  GENERAL  THEORY  OF  ELECTRO¬ 
MAGNETIC  MODES  OF  A  CYLINDER 


The  solution  to  the  electromagnetic  modes  of  a  cylinder 
has  been  known  for  quite  some  time26,27  and  is  briefly  de¬ 
rived  here. 

We  consider  a  cylinder  of  radius  R  of  dimensionless  elec¬ 
tric  permittivity  e2,  centered  along  the  z  axis  and  surrounded 
by  a  second  dielectric  medium  et.  For  nonmagnetic  media, 
the  electric  and  magnetic  fields  in  frequency  space  satisfy  the 
wave  equation 


V 


E(r) 

H(r) 


E(r) 

H(r) 


=  0. 


(Al) 


We  have  demonstrated  that  the  subwavelength  confine¬ 
ment  of  guided  plasmon  modes  on  conducting  nanostructures 
leads  to  strong  coupling  between  these  modes  and  nearby 
emitters  in  the  optical  domain.  This  strong  coupling  leads  to 
large  effective  Purcell  factors  for  emission  into  the  plasmon 
modes,  which  are  limited  only  by  heating  losses  in  the  con¬ 
ductor.  While  losses  prevent  the  plasmon  modes  from  being 
useful  photonic  carriers  of  information,  we  have  shown  that 
they  can  be  efficiently  outcoupled,  e.g,  to  a  dielectric  wave¬ 
guide.  We  estimate  that  single-photon  generation  efficiencies 
exceeding  95%  are  possible  in  such  a  tiered  system.  Finally, 
we  have  analyzed  the  effects  of  plasmon  scattering  due  to 
moderate  surface  roughness  on  these  nanostructures. 

Rapid  advances  in  recent  years  in  fabrication  techniques 
for  nano  wires, 49,50  nanotips,51  and  subwavelength  dielectric 
waveguides52,53  put  such  a  system  in  experimental  reach. 
Quantum  dots  or  single  color  centers  might  serve  as  physical 
realizations  of  solid-state  emitters,  which  could  be  used  to 
achieve  strong-coupling  cavity  QED  and  quantum  informa¬ 
tion  devices  on  a  chip  at  optical  frequencies.  It  is  also  inter¬ 
esting  to  consider  real,  individual  atoms  interacting  with 
nanowires  and  the  challenges  associated  with  constructing 


The  solutions  to  Eq.  (Al)  are  perhaps  most  easily  derived  by 
first  finding  scalar  solutions  of  the  equation  and  then  con¬ 
structing  vector  solutions.  Working  in  cylindrical  coordi¬ 
nates,  scalar  solutions  of  Eq.  (Al)  satisfying  the  necessary 
boundary  conditions  take  the  form  if/1^Hm(ki_Lp)eim^+lkwz 
and  02oc^m(&2±pVm<^HA:||Z  outside  and  inside  the  cylinder,  re¬ 
spectively.  Here,  Jm  and  Hm  are  Bessel  functions  and  Hankel 
functions  of  the  first  kind,  respectively,  while  kil  =  ^kf-kfl 
and  kt=  o>V^/ c .  Jm  is  well  behaved  at  p= 0,  while  Hm(x ) 
~  elx  for  large  v  satisfies  outgoing-wave  conditions.  It  is  easy 
to  verify  that  two  independent  vector  solutions  to  Eq.  (Al) 
are  given  by  vp^VX  (£ 0-)  and  w*=  \  V  X \t.  The  curl  rela¬ 
tions  of  Maxwell’s  equations  then  imply  that  E  and  H  must 
take  the  form 

E;(r)  =  a(v((r)  +  />,w,(r) ,  (A2) 

H,(r)  =  -  —k,[a,vf,(r)  +  /?,v,(rj] ,  (A3) 

(O/X0 

where  at  and  bt  are  constant  coefficients.  Expanding  out 
these  expressions  in  detail, 
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E,(r)  = ' 
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(A4) 


where  FUm(x)=Hm(x)  and  F2m(x)=Jm(x). 

Up  to  this  point  at  and  bt  are  arbitrary  coefficients,  whose 
relationship  becomes  fixed  by  imposing  boundary  conditions 
between  the  cylinder  and  surrounding  dielectric.  Requiring 
that  the  tangential  components  E Ev  H and  Hz  of  the 
fields  be  continuous  at  the  boundary  results  in  a  linear  sys¬ 
tem  of  four  equations,  which  we  write  in  abbreviated  matrix 
form  as  M{ax  a2  bx  b2)T=  0.55  A  nontrivial  solution  for  the 
fields  requires  that  detM=0,  which  after  some  work  simpli¬ 
fies  to  the  mode  equation  given  in  Eq.  (1). 

One  special  case  of  interest  is  that  of  a  TM  mode  with  no 
winding  (m  =  0).  The  component  of  H  along  £  by  definition 
vanishes,  which  implies  that  the  coefficients  at  in  Eq.  (A4) 
vanish.  The  condition  detM=0  is  then  significantly  easier  to 
evaluate  in  this  situation.  In  particular,  at= 0  implies  that  the 
field  components  E ^  and  Hz  vanish,  and  continuity  of  the 
remaining  tangential  field  components  Ez  and  H ^  at  the 
boundary  requires  that 


/  t l 

7 THo(ki±R) 

Kl 

—k1±H'0(kl±R) 

\  WyU,0 


-  -jrUhxR) 


bx 

^2 


(A5) 


Setting  the  determinant  of  the  above  matrix  equal  to  zero 
immediately  yields  the  mode  equation  of  Eq.  (2),  and  it  is 
also  immediately  seen  that  the  ratio  of  the  coefficients  bX  2 
must  be  given  by  Eq.  (4). 


APPENDIX  B:  DERIVATION  OF  CUTOFF 
FOR  HIGHER-ORDER  MODES 

In  this  section,  we  show  that,  to  a  very  good  approxima¬ 
tion,  a  nano  wire  essentially  supports  a  single,  fundamental 
m  =  0  plasmon  mode.  In  particular,  for  all  higher-order  plas- 
mon  modes  |m|^2,  a  cutoff  wire  size  Rcutoff  exists  below 
which  such  modes  cannot  exist,  while  the  |m|  =  l  plasmon 
modes  exhibit  an  exponential  growth  in  their  mode  volumes 
as  R— » 0.  For  simplicity,  we  will  assume  in  this  section  that 


we  are  dealing  with  a  lossless  system  (Im  62=0). 


1.  Behavior  of  \m\^2  modes 

We  are  interested  here  in  the  behavior  of  the  |m|^2 
modes  near  cutoff,  which  is  characterized  by  a  small  devia¬ 
tion  of  the  plasmon  wave  vector  k\\  from  V^w/c  (see  Fig.  2). 
To  simplify  algebra  in  the  derivation  of  Rcutoff,  from  this 
point  forward  we  make  the  mode  Eq.  (1)  dimensionless  by 
setting  (x)/c=  1,  and  we  will  assume  that  m  is  positive  (the 
case  where  m  is  negative  follows  this  derivation  with  a  few 
minor  modifications).  Furthermore,  it  is  useful  to  define  a 
small  quantity  S=k\\-Je[,  where  we  specifically  consider  the 
positive  k\\  solution.  On  physical  grounds,  any_mode  with 
positive  k\\  must  satisfy  8^ 0,  because  if  k\\<iex,  the  fields 
outside  the  wire  would  be  radiative  in  nature  and  implies  that 
the  system  is  continually  radiating  energy  out  to  infinity 
without  a  source.  It  follows  that  any  value  of  R  where  £=0 
becomes  a  solution  to  Eq.  (1)  for  some  m  then  corresponds  to 
a  critical  point  in  behavior,  and  specifically  is  a  cutoff  be¬ 
yond  which  modes  cease  to  exist  for  that  given  m.  To  find 
this  R=Rcut0 ff,  we  expand  the  two  sides  of  Eq.  (1)  in  8.  Both 
sides  have  contributions  to  these  expansions  that  are  diver¬ 
gent  at  £=0  (terms  that  behave  like  8~n ,  where  n>  0),  and 
we  will  show  that,  for  m  ^  2,  there  exists  one  value  of  R  that 
equates  these  two  divergent  contributions;  i.e.,  £=0  satisfies 
the  mode  equation  at  this  particular  value  R=Rcntoff- 

It  is  straightforward  to  show  that  the  divergent  contribu¬ 
tion  to  the  expansion  of  the  left-hand  side  (LHS)  of  Eq.  (1)  is 
given  by 


LHS  = 


m 2  m2( 361  +  62) 

4 R2  82  4 R2  yje1  ( ex  -  62)  8 


+  e>(<5°). 


(Bl) 


To  expand  the  right-hand  side,  we  first  note  that  the  quan- 
tity  (1  /k2±)J'm(k2±R)/Jm(k2±R)  =  (1  /  ^2-eM^2-e,R) 
+  (D(8l)  is  well  behaved  near  £=0.  Here,  we  have  defined 

Jm(x)=J'm(x)/Jm(x).  Then,  using  the  identity 


Elm {ix)  rn+ 1  (-^ )  > 

777 


(B2) 


where  Km(x )  is  a  modified  Bessel  function  of  the  second 
kind,  and  the  expansions 


Km{x)  = 


(m  -  1)! 


(m  -  2) ! 


^  \  m—2 
X) 


+  0(x4~m) 


(m  S'  2), 

Sm 


2mem 


+  0{$a) 


(B3) 

(B4a) 


=  iKl±,  (B4b) 

it  is  tedious  but  straightforward  to  expand  the  expression 
(1  /k1_L)[H'm(kilR)/Hm(kllR)]  as  well.  Performing  these  ex¬ 
pansions  and  simplifying,  one  finds  that 
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RHS  = 


4  R2S2 


m 


4  Wei  <5 


mVe i” 

2(m-  1)<5 


/m(61  +  g2)-/m(Vg 2~gl^) 

2RVe!(e!  -  £2)<S 


■0(<5°).  (B5) 


Comparing  Eqs.  (Bl)  and  (B5),  we  see  that  S=0  is  a  solution 
provided  that  these  terms  are  equal  to  0(S~l),  i.e., 


me i  +  e2  _  Ret  I(ei  +  e2)Jm(Ve2 -  6tR) 
R  e2  -  ex  m  -  1  Ve,  -  e2 


(B6) 


The  solution  /? = /?cutoff  to  Eq.  (B6)  gives  the  cutoff  wire  size, 
below  which  the  mode  m  cannot  exist.  In  the  regime  of  in¬ 
terest  (ei>0,  62 < 0,  61  +  62<0),  the  first  and  second  terms 
are  positive  while  the  third  term  is  a  negative  function  (for 
R > 0)  that  behaves  like  -HR  for  small  R  and  approaches  a 
constant  for  large  R.  It  can  be  seen  then  that  a  solution  exists 
for  any  m  ^  2,  which  establishes  that  these  modes  are  indeed 
cut  off  in  the  nanowire  limit. 


2.  Behavior  of  |m|  =  l  mode 

For  simplicity,  we  will  assume  that  m- 1,  as  the  case  of 
m  =  - 1  follows  this  derivation  closely.  The  case  of  m=  1  must 
be  studied  separately  because  the  expansion  of  Km(x)  given 
in  Eq.  (B3)  only  holds  for  2.  The  different  asymptotic 
scaling  of  Kx(x)  leads  to  a  unique  behavior  of  the  m- 1  mode 
in  the  nanowire  limit.  In  particular,  we  will  show  that  this 
mode  does  not  strictly  have  a  cutoff  size,  but  that  k\\  — ► 
exponentially  in  the  limit  R^  0.  In  turn,  the  magnitude  of 
klx  becomes  exponentially  small,  which  corresponds  to  an 
exponential  growth  in  the  spatial  extent  or  mode  volume. 

Again  defining  S=k we  are  interested  in  finding  an 
approximate  solution  to  Eq.  (1)  in  the  limit  of  small  S.  We 
proceed  by  expanding  both  sides  of  the  equation  as  a  series 
in  the  small  parameter.  The  expression  for  the  left-hand  side 
given  by  Eq.  (Bl)  remains  valid  for  m=  1.  For  the  right-hand 
side,  we  anticipate  that  both  the  quantities  kl±R  and  kl±  will 
be  small  as  R— >0  (these  assumptions  can  be  checked  for 
consistency  at  the  end  of  the  calculation),  and  we  thus  ex¬ 
pand  around  k1±R  =  0  the  term 

i  g;(*1±*)_  i  (B7a) 

k\LHi(hiR)  iKxLHx{iKlLR) 


1 


iKlA><lLR 
+  0{k\lR2) 


l  .(  ,  Kl±R 

-i\  7+ log  — — 


Kl±R 


(B7b) 


where  0.577  is  Euler’s  constant.  Here,  we  have  used  Eq. 
(B2)  to  convert  Hm(ix)  to  Km(x)  and  the  expansion 


1  HlikxJi)  1  /  .  k1±R\ 

ktlHt(kt±R)  iKll\K1±R  lKl±  °8  2  /’ 

(B9) 


Furthermore,  having  assumed  that  k1±  (and  by  extension, 
/c1±)  is  a  small  quantity,  we  can  now  expand  the  expression 
above  in  terms  of  S  using  Eqs.  (B4a)  and  (B4b).  Making  this 
substitution,  and  after  a  bit  of  algebra,  one  finds  that  the 
expansion  of  the  right-hand  side  (RHS)  of  Eq.  (1)  is  given  by 


„TTO  _  1  el  +  3e2-2R2el(e1~  e2)log(SR2'J7l/2) 

RHS  —  0  0  „  i — 

4 R  S  4R2Vw(6i  -  62)^ 

(B10) 

Finally,  equating  the  left-  and  right-hand  sides  to  0(8~l) 
gives  the  solution 


2(6i  + 

^6l(62  _  6l)  _ 


(Bll) 


It  follows  that  in  the  nanowire  limit, 

KlL  =  {k]-ex)m  (B12a) 


=  (2<5Ve^)1/2  (B12b) 


2 

~  —  exp 
R  F 


61  +  62 

^6l(62  _  6l)  _ 


(B12c) 


Equations  (Bll)  and  (B12c)  indicate  that  the  m=  1  plas- 
mon  mode  does  not  have  a  cutoff  in  the  nanowire  limit,  but 
instead  that  its  longitudinal  wave  vector  approaches  ex¬ 
ponentially,  with  a  corresponding  exponential  increase  in  its 
transverse  extent  (~l//c1±)  and  mode  volume.  It  is  therefore 
well  justified  to  say  that  this  mode  is  effectively  cut  off,  as 
the  coupling  strength  to  this  mode  becomes  strongly  sup¬ 
pressed  as  R^0. 


APPENDIX  C:  RADIATIVE  AND  N ONRADI ATI VE 
DECAY  RATES  NEAR  A  NANOTIP 


Here,  we  derive  more  carefully  the  expressions  given  in 
Eqs.  (41)  and  (42)  for  the  radiative  and  nonradiative  sponta¬ 
neous  emission  rates  near  a  nanotip. 

To  calculate  the  radiative  rate,  we  consider  our  expression 
for  <Er  in  Eq.  (36)  in  the  far- field  (large  v )  limit,  where  the 
Km(qv)  terms  in  Or  decay  exponentially  with  v.  Because  of 
this  exponential  dependence  at  large  v ,  to  good  approxima¬ 
tion,  it  suffices  to  expand  the  terms  am{q),  Jm(qu '),  and 
Km(qv')  around  q  =  0.  The  only  nontrivial  expansions  occur 
for  the  terms  am(q)  and  are  given  by 

«o(<?)  =  +  (Cl) 


log  2 
2 


1  log  V 
4  +  2 


I  x  +  0{x2). 


(B8) 


aM)  =  ~~ — —  q2vl+0(q4),  (C2) 

2  £\  +  e2 


We  now  assume  that  kx±R  is  small  enough  that  y 
<^|log  k1±R|,  such  that 


and  so  on.  These  expansions  allow  for  exact  evaluations  of 
the  integral.  It  can  be  verified  that  the  dipole  contributions  to 
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Or  originate  from  the  m  =  0,l  terms  in  the  sum,  which  are 
readily  found  to  be 


<J>(m=0)(r,r') 


1 _  2(  _  e2\4(t;2-M2) 

477e0e1l,°\  ej(u2  +  v2)3 


In  +  <5<E(r), 
v 


(C3) 


^m=1)(r,r') 


1 


47T6n6 


cos(</>-  (j)') 


61-62  ^0^ 


ofci 


6i  +  62  l/  (U2  +  U2)3 

(C4) 


Here,  is  a  complicated  function,  but,  most  importantly, 
contains  no  dependence  on  r'.  Recalling  that  the  pseudopo¬ 
tentials  derived  above  correspond  to  a  point  charge  source, 
we  can  immediately  obtain  the  potentials  due  to  a  dipole 
\)0e~lO)t  at  r'  by  applying  the  operator  (p0-V')  to  these  ex¬ 
pressions.  In  parabolic  coordinates,  the  gradient  operator  is 
given  by 


V  = 


Vw2  +  v2 


„  .  1  ,  d 

u —  +  v —  +  — (b — 7, 
du  dv )  uv  dcf) 


(C5) 


and  for  a  dipole  located  on  the  z  axis  (u'  =  0),  we  find  that 


<D 


(m=0)  __ 
dip,r 


1 

47T£()€\ 


^2 \  vl  (#  ~  r)(Pp  '  v) 
ejv'2  r3 


(C6) 


,  1  gi  -  g2  vl  [po  -  g(Po  •  g)]  •  r 

dlp’r  4tt€0£i  6[  +  e2  v' 2  r3 

From  these  expressions,  one  can  immediately  identify  the 
induced  dipole  moments  in  the  nanotip, 

^=-^0772(1  -  7^)  (p<>ii£)>  (C8) 


£p  =  up0— - -  (Po  i-  £) ,  (C9) 

61  +  62  U 

and  arrive  at  the  radiative  decay  rates  given  in  Eq.  (41). 

The  leading  term  for  the  nonradiative  decay  rate  is  found 
by  calculating  the  divergence  in  the  reflected  field  Er(r'  ,r') 
as  vf  — >v0.  The  reflected  field  Er=-V(p0- V')<Fr  is,  in  gen¬ 
eral,  difficult  to  evaluate,  but  simplifies  considerably  for  a 
dipole  located  on-axis  (u'=  0)  due  to  the  presence  of  the 
Jm{qu’)  term  in  Or,  given  in  Eq.  (36).  The  operation  V' 
causes  terms  like  Jm( 0)  and  J'm( 0)  to  appear  in  Er,  which  are 
nonzero  only  when  m  =  0  and  m=  1,  respectively.  This  imme¬ 
diately  leads  to  the  expressions 

2  r  00  3 

Po'Er(r',r')  =  -  Po  dq^ax{q)K\(qv')  (p0-L£), 
4T7e0eiJ0  vz 

2  roc  3 

p0-Er(r',r')  =  -  P°  dq-^a0(q)K2(qv')  (p0ll£), 

27j-e0£iJo  v 

(CIO) 

which  were  given  in  Eq.  (40).  Examining  further  the  solu¬ 
tions  to  a01,  it  can  easily  be  shown  that  their  asymptotic 


expansions  in  the  limit  qv0>  1  take  the  form 

“0,1  (?)  “  - - -e2qv°  (qv 0>1).  (Cll) 

77  61  +  62 

At  the  same  time,  in  the  limit  qv'  >  1 ,  the  behavior  of  K\  is 
given  by  K\(qv')~{7r/2qv')e~2qv  ,  and  thus,  as  v'  ^v0  the 
integrands  of  Eq.  (CIO)  exhibit  very  long  tails  due  to  the 
presence  of  terms  ~e~2q<<v  ~Vo)  at  large  q.  The  tail  is  the  origin 
of  the  divergence  that  we  expect  on  physical  grounds.  Using 
these  expansions  as  well  as  the  fact  that  the  decay  rate  is 
proportional  to  Im(p0-E),  the  integrals  can  be  evaluated  ex¬ 
actly  and  yield  the  nonradiative  decay  rates  given  in  Eq. 
(42). 


APPENDIX  D:  BOUNDARY  ELEMENT  METHOD 

Our  numerical  implementation  of  the  boundary  element 
method  (BEM)  closely  follows  the  method  derived  in  Ref. 
40.  Here,  we  briefly  outline  the  main  ideas  of  BEM  while 
referring  the  reader  to  Ref.  40  for  more  details,  and  we  dis¬ 
cuss  the  key  elements  of  our  implementation. 

We  assume  that  our  system  contains  a  set  of  known,  time- 
harmonic  source  charges  and  currents  pext  Jext  in  the  pres¬ 
ence  of  some  scattering  dielectric  body  whose  surface  is  de¬ 
noted  S  (although  we  discuss  one  body  here,  BEM  is  easily 
generalizable  to  treat  multiple  scatterers).  In  the  case  of  in¬ 
terest,  S  represents  the  surface  of  a  metallic  nanotip,  while 
the  external  source  corresponds  to  an  oscillating  point  dipole 
\)0e~l(Ot  at  some  location  r\  For  simplicity,  we  also  assume 
that  we  are  working  with  nonmagnetic  media,  and  we  denote 
by  €j  (7=  1,2)  the  dimensionless  electric  permittivities  out¬ 
side  and  inside  S ,  respectively.  The  underlying  principle  be¬ 
hind  BEM  is  that  the  scalar  and  vector  potentials  and 
A;(r)  in  each  region  can  be  written  (in  the  Lorentz  gauge)  in 
the  form 


w*) 


=  —  \  dr’Gj(r  -  r')pext(r')  +  -  - —  f  < 

4TT€0ej  J  4ire0£jJs 


XG,(r-s)cr,(s), 


ds 


(Dl) 


A;(r)  =  ~~  I  dr’ G4r  -  r')jext(r')  +  7-^  I  dsG,(r-s)h,(s), 


477 


mo  r 

477-Jc 


(D2) 


ikjr 

Gj(r)  = - ,  (D3) 

r 

where  Gj  is  the  Green’s  function  in  a  medium  of  uniform  ep 
and  kj=4ej((x)/c).  Physically,  the  equations  above  state  that 
the  fields  in  region  j  can  be  described  as  a  result  of  the 
combination  of  the  external  sources  and  some  effective  sur¬ 
face  charge  and  current  distributions  ctj  and  h;  on  S.  In  gen¬ 
eral,  these  effective  distributions  do  not  have  physical  sig¬ 
nificance;  for  example,  they  do  not  correspond  to  actual 
charges  and  currents,  and  the  distributions  in  region  1  and 
region  2  are  not  necessarily  equal  [e.g.,  ax( s)  A  (r2(s)].  The 
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values  of  dj  and  hj  are  not  known  initially,  but  a  set  of  linear 
integral  equations  for  these  distributions  results  from  enforc¬ 
ing  various  boundary  conditions  for  the  scalar  and  vector 
potentials  at  S.  In  particular,  A,  Dx,  and  Hy  must  be 
continuous  at  the  boundary.  To  calculate  the  distributions  nu¬ 
merically,  if  the  boundary  S  is  finite,  one  can  mesh  up  the 
surface  into  a  finite  number  of  grid  points.  Assuming  that  crj 
and  hj  are  constant  over  each  grid  point,  the  linear  integral 
equations  become  a  set  of  linear  equations  in  the  values  of  crj 
and  hj  that  can  be  solved  straightforwardly.  Once  these  dis¬ 
tributions  are  known,  the  potentials  and  then  the  fields  E  and 
H  can  be  calculated. 

In  our  problem  of  interest,  we  assume  that  the  dipole  is 
located  on  the  z  axis  and  oriented  along  £,  while  the  nanotip 
is  described  by  a  paraboloid  of  revolution  around  the  z  axis. 
Due  to  the  axial  symmetry  of  the  system,  BEM  simulations 
are  advantageous  because  one  only  needs  to  calculate  the 
unknown  distributions  along  one  dimension  instead  of  over 
the  entire  two-dimensional  surface  S.  At  the  same  time,  the 
source  is  a  dipole  oscillating  at  constant  frequency,  and  thus, 
the  external  charges  and  currents  are  calculated  quite  easily. 
In  BEM  (at  least  in  the  current  formulation),  it  is  necessary 
that  the  nanotip  surface  S  be  finite,  and  we  implement  this 
numerically  by  tapering  and  rounding  off  the  nanotip  far 
from  the  region  of  interest.  In  general,  any  termination  can 
result  in  some  back  reflection  of  the  guided  plasmon,  and  this 
results  in  some  small  oscillations  of  the  fields  due  to  inter¬ 
ference  with  the  forward  propagating  plasmon,  as  barely 
seen,  e.g,  in  Fig.  5.  In  our  simulations,  the  reflected  ampli¬ 
tude  is  kept  to  within  a  few  percent.  Very  fine  meshes  were 
used  to  ensure  accuracy;  in  most  of  our  simulations,  for  ex¬ 
ample,  the  spacing  between  points  in  the  regions  of  constant 
R  was  approximately  V/400. 

APPENDIX  E:  DERIVATION  OF  COUPLED-MODE 
EQUATIONS 

In  this  section,  we  derive  the  equations  of  evolution  for 
two  electromagnetically  coupled  systems  based  on  Lorentz 
reciprocity. 

First,  we  derive  the  Lorentz  reciprocity  equation  gener¬ 
ally.  Assuming  nonmagnetic  media,  suppose  that 
{E^r^H^r^e^r)}  and  {E2(r)  ,H2(r) ,  62(r)}  separately  sat¬ 
isfy  Maxwell’s  equations.  At  this  point,  the  systems  1  and  2 
and  their  field  solutions  are  not  necessarily  related  to  each 
other  at  all.  In  the  following,  we  assume  that  all  fields 
E(r,t)=E(r)e~l0)t,H(r,t)=H(r)e~lb)t  have  harmonic  time  de¬ 
pendence.  Using  the  vector  identity 

V  •  (a  X  b)  =  b  •  (V  X  a)  -  a  •  (V  X  b)  (El) 
and  the  curl  relations  of  Maxwell’s  equations,  we  can  write 

V  •  (Ej  X  H*)  =  H2  •  (V  X  Ej)  -Er(VX  H*) 

=  H;  •  (;w/z0H|)  -  E,  •  (icoe0e2E*2),  (E2) 

and  similarly, 
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V  •  (E2  X  Hj)  =  Hr(VX  E*2)  -  E*2  •  (V  X  Hj) 

=  H!  •  (-  icop0 H2)  —  E2  •  (—  iu>e oeiEi) . 

(E3) 

Adding  up  Eqs.  (E2)  and  (E3)  yields  the  equation  for  Lorentz 
reciprocity, 

V  •  (Ej  X  h;  +  E*2  X  HO  =  i(O€0Ex  ■  E^r)  -  e2(r)]. 

(E4) 

We  now  derive  coupled-mode  equations  for  two 
waveguides  based  on  the  Lorentz  reciprocity  equation  above. 
This  derivation  closely  follows  that  of  Ref.  46.  We  empha¬ 
size  that  the  nature  of  the  waveguides  can  be  quite  general; 
e.g.,  they  can  be  any  type  of  normal  dielectric  or  plasmon 
waveguide.  We  let  the  indices  jm,v=a,b  refer  to  the  system 
consisting  of  waveguide  a  without  the  presence  of  system  b, 
and  b  the  system  consisting  of  waveguide  b  without  the  pres¬ 
ence  of  a.  We  also  assume  that  the  surrounding  dielectrics 
for  systems  a  and  b  are  the  same,  i.e.,  ea(r=^)  =  €b(r=^), 
and  that  the  waveguides  are  copropagating  along  the  z  direc¬ 
tion.  It  is  assumed  that  the  total  electric  field  for  the  system 
consisting  of  waveguides  a  and  b  together  can  be  written  as 

Er(r)=  2  C„(z)E,(r),  (E5) 

v=a,b 

with  a  similar  expression  for  H.  That  is,  we  assume  that  the 
total  field  can  be  written  as  a  linear  superposition  of  the 
unperturbed  modes  of  systems  a  and  b.  For  the  case  where 
systems  a  and  b  each  have  one  allowed  mode,  the  index  v 
refers  to  these  unperturbed  modes.  In  general,  when  a  and  b 
have  Na  b  allowed  unperturbed  modes,  v  is  understood  to  be 
an  index  that  covers  all  of  these  modes.  We  can  derive  exact 
equations  of  motion  for  Cv(z)  by  using  Eq.  (E4).  Specifically, 
we  will  let  the  index  1  =  T  in  Eq.  (E4)  refer  to  the  total  fields 
Er(r)  and  Hr(r)  and  the  dielectric  profile  of  the  combined 
system  6r(r),  while  we  will  let  the  index  2  =  /jl  refer  to  any 
one  of  the  allowed,  unperturbed  modes  of  systems  a  and  b. 
Substituting  this  into  Eq.  (E4)  yields 

V  •  (Er  X  +  E^  X  Hr)  =  ico€0ET  •  E^[6r(r)  -  ^(r)] 

(E6) 

or 

2  V  •  [C„(z)E„  X  h;  +  C,(z)e;  X  HJ 

v=a,b 

=  ic*e0  2  Cv{z)Ev  •  E;Mr)  -  e*(r)] .  (E7) 

v=a,b 

Applying  the  divergence  theorem  to  this  result  (and  assum¬ 
ing  that  the  surface  terms  vanish  at  infinity)  gives 

7-2  I  d2p[Cv{z)Ev  X  h;  +  C„(z)E;  x  H  J  •  z 

v=a,b  2 

=  ioco  2  C„(z)  f  d2pEv  ■  E^[er(r)  -  e*  (r)].  (E8) 

v=a,b  J 

The  left-hand  side  can  be  further  simplified, 
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LHS 


=  j 2  f 

°Zv=a,b  J 


d2p[C„(z) E„  X  H*  +  C„(z)E*  XHj-f 


(E9a) 


_  V*  ClC  y 

—  ^  rjy 

v=a,b  ^ 


+  Cj/z) 


J  d2P{  e,xh;  +  e;xh,)-£ 

cPp(E„  XH’+  E*  X  H„)  •  £ 


(E9b) 


r/C, 


=  2  * 

v=ci,b 


/■ 


d2p(Ev  X  +  E^  X  H„)  •  £ 

+  C,/z)(  ;oje0  J  d2pEv  •  E*[ev(r)  -  e*(r)]  j ,  (E9c) 


E  total _ 

2  “ 


^ife-L  ,  w 


Mofe_LPVi||Z 

n2 

r°°  ih  /i 

+  pl  l^221-/o(/i2-LP)g'/‘"Zg(^l|) 

7  _c 


'  —00 

,2 


^fyb2J0(k2±p)e,k^ 
k2 


’[ 


h 


+  p\  dh]l-^rJ0(h21p)e,tl"zB(hl) 


z, 


^ytotal _ 

1,0 


//ci  Lb\H'(){k\ Lp)e,kvz 

roc 

+  p  I  dh\\ihl_LHl)(hi1p)elhuzA(h^ 

J  —00 


where  we  have  applied  the  divergence  theorem  on  Eq.  (E4) 
to  get  the  last  line. 

Substituting  Eq.  (E9c)  back  into  Eq.  (E8)  yields  a  set  of 
Na+Nb  coupled,  first-order  differential  equations: 

X  ^pPVM{z)  =  -  iwe0X  Cviz)Km[z), 

v  v 

(E10) 

PvJz)  =  J  d2p( E„  XH’+E^x  Hv)  •  f , 

(Ell) 

KVJU(z)  =  J  d2pEr  ■  E*[e„(r)  -  e7(r)]. 

(El  2) 

We  emphasize  that  these  coupled-mode  equations 
within  the  ansatz  given  by  Eq.  (E5). 

are  exact 

APPENDIX  F:  RADIATIVE  SCATTERING 

Suppose  that  in  the  presence  of  roughness,  the  first-order 
scattered  field  Ez  is  given  by  Eq.  (59).  Using  the  expressions 
derived  in  Eq.  (3)  for  the  unperturbed,  incident  plasmon 
field,  and  letting  k\\  denote  the  unperturbed  plasmon  wave 
vector,  the  total  (incident  plus  scattered)  fields  to  first  order 
in  p  are  given  by 


g total  _ 


l^^btH[>(k\1p)e,k*z 


+  p  |  dh{fifpH'0(hllP)e^A(h^ 
k\ 


~72'blH0(kllp)e'k"z 


+  p  dh\\—yH0{hlLp)el^zA{hy) 


z. 


H 


rtotal 

2,0  : 


1 

«yU0 


ik2±J'0{k21p)b2e,k"z 


+  p  dhllih21Jo(h21p)elh"zB(hll) 


(FI) 


The  boundary  condition  equations  in  Eq.  (60)  can  be 
solved  by  plugging  in  the  fields  above,  carefully  expanding 
the  equations  as  a  power  series  in  p,  and  then  solving  for 
each  order  of  p,  utilizing  the  expansions 

t  =  z  +  p~fp  +  0(p2), 


FiJki±p0)  =  Flm(kuR)  +  plkl±F’Jkl±R)  +  Oip2), 
Fi,m(hi±p o)  =  FlJhilR)+p£hilF’iJ.hllR)  +  Oip2), 

(F2) 


where  Flm(x)=Hm(x)  and  F2m(x)=Jm(x).  The  resulting 
O(p0)  equations  are  trivially  satisfied  by  the  plasmon  fields 
of  a  smooth  nanowire,  while  the  0(p)  equations  are  found  to 
be 


f 


dhw 


^k//0(/j1±i?)A(^ii)-  — 


J0(h2±R)B(h^ 


oihf- 


X2± 


b2az)Jo(k21R )  -  -±b^iz)H'0iklLR) 


ik\\k2Li  di  ihkii,  d£ 


b2fJo(k2±R)~- 

dz 


b^H'0iklLR) 

dz 


,lk¥- 


dhlhlLH'0{hlLR)A{h^  -  h2LJ'Q(h21_R)B{h^eih*z 

oo 

=  [k221b2az)fo(k21R)  -  klMzmk^Rfte^. 

(F3) 


Here,  we  assume  that  the  metal  inherently  has  no  losses,  i.e., 
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Im  62— 0,  such  that  k\\  is  purely  real.  Then,  by  plugging  in  the 
Fourier  transform  of  £(z)  given  in  Eq.  (61),  the  equations 
above  become  purely  algebraic.  It  is  tedious  but  straightfor¬ 
ward  to  show  that  the  solutions  are  given  by  Eq.  (62),  with 
the  coefficients  /(fy|)  and  g(h^  defined  via 


m = 


g(h  n)  = 


h2±NJ0(h2±R)  -  ^MQi^J^R) 

H0{hlLR)J'02LhlLk\  -  H^h1±R)J0(h2±R)h2±k2  ’ 

hl±NH0(hlxR)  -  kjMjh^H^h^R) 
H0(hl±R)Jfi±hl±%  -  H'0(h1±R)J0(h21R)h2±lZ  ’ 


hki±(h\\~k\\) 


hK(kixR) 


k\\k2±(h\\-k^ 


ViM). 


N  =  b1k\±HlQ(k1±R)  -  b2kljl(k2LR) .  (F4) 

To  evaluate  the  expression  for  rradrough  in  Eq.  (68),  we 
first  normalize  the  fields  such  that  the  denominator  takes  the 
value  h(o/4L.  For  concreteness,  we  assume  that  the  conduct¬ 
ing  nanowire  exhibits  a  Drude-like  behavior  and  that  we  are 
operating  well  below  the  plasma  frequency,  such  that  e2((o) 
=  1  -  up  co2  ~  -(Op/  co2 .  In  this  case,  the  dispersive  term  in  the 
denominator  is  positive  and  given  by  ^[co62(co)]~  |e2(a>)|. 
In  the  nanowire  limit,  normalization  of  the  integral 
is  straightforward  and  yields  a  coefficient  bx 


~  ^hcok^^R2 /  €qVC4L.  Here,  V  is  a  dimensionless  parameter 
given  by 

'(C) 

-T—  dxx[I](x)  +  /q(x)] 

U(C) Jo 


-  8e?  (  K2AC) 

V--  11  1 


nC2\\e2\ll 


+  —J  dxx{K\(x)  +  K^(x)\ ) . 


(F5) 


Then,  using  the  relationships  k\\~C/R  and  ki±~iC/R ,  one 
can  calculate  the  leading  terms  of  as  R  — ►  0, 


h22bl{CIR)2<f> 

-  2ik\h2J irhi  lR  ' 


(F6) 


where  we  have  defined  (p=  (b2l b AJ'^iiC) - H'^dC) .  In  the 
equation  above,  we  have  explicitly  given  the  leading  terms 
of  the  numerator  and  denominator  of/(/zy).  The  ratio  Z?1/Z?2  is 
given  in  Eq.  (4),  which  in  the  nano  wire  limit  results  in  the 
simplification  of  (j)  given  in  Eq.  (70). 


APPENDIX  G:  N  ONR  ADI  ATI  VE  SCATTERING 


The  elements  of  the  matrices  Mt  and  vectors  yt  appearing 
in  the  matrix  integral  Eq.  (79)  in  the  presence  of  surface 
roughness  are  given  by 


/  KjhR)  - IJhR )  \ 
\heiK'm(hR)  -  heJ'JhR)  ) 


(Gl) 


J 


Ml(h,q)=  I 


ciK'JqR) 
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J{h-q)z 
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(G2) 
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(G3) 


Mf(h,q,q')  =  ex 


a3  ~  1 

^KZ(qR)  +  q2(h  -q-  q’)K'm(qR)  +  -qq'Qi  -q-  q')K'm{qR ) 


(G4) 


M22(h,q,q')  =  -  e2 


a3~  1 
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v0(A) =  - 
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/ 


\2(h,q,q')  =  - 


\ 


\2 


/"'(^)  +  <r(^ -<?- q')I'm(qR)  +  ~qq'(h  -q-  q')I'm(qR) 


Kjqp')ei(h^' . 


/ 


(08) 
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